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CHAPTER  1 


INTRODUCTION 


1.1  General 

Heavy  metals  in  our  aquatic  environments,  their 
transport  processes,  env ironmental  distributions  and 
biological  effects  have  been  the  focus  of  increasing 
concern.  Direct  toxicity  to  aquatic  organisms  and  in¬ 
direct  toxicity  to  higher  organisms  and  man  are  the 
center  of  this  concern.  Although  heavy  metals  are  natural 
constituents  of  aquatic  environments,  additional  quan¬ 
tities  introduced  by  man's  presence  may  result  in  an 
altered  chemical  composition  of  the  aqueous  and  sediment 
phases,  often  with  detrimental  impacts  on  the  environ¬ 
ment.  in  an  effort  to  reduce  these  impacts,  environmental 
transport/fate  models  of  heavy  metals  in  aquatic  ecosys¬ 
tems  are  being  developed. 

The  purpose  of  this  study  was  to  evaluate  the 
Metal  Exposure  Assessment  Modeling  System  (MEXAMS).  This 
program  is  the  result  of  linking  MINTED,  a  geochemical 
model,  and  the  Aquatic  Exposure  Assessment  Model,  EXAMS. 
MEXAMS  was  developed  for  the  Env i ronmental  Protection 
Agency  (EPAj  (Felmy  et  a!.,  1983)  to  allow  assessment  of 


the  impacts  of  priority  pollutant  metals  (As,  Cd,  Cu,  Pb, 
Ni,  Ag  and  Zn).  The  evaluation  consisted  of  calibrating, 
testing,  and  assessing  the  model  response  for  actual  site 
application  on  the  Naugatuck  River  located  in  the 
Housatonic  River  Basin  in  western  Connecticut. 

1.2  Background  Information 

The  Naugatuck  River  is  located  in  western  Connec¬ 
ticut  and  has  a  long  history  of  industrialization  along 
its  length.  Hydropower  generation  was  responsible  for  the 
initial  development  of  the  valley,  with  light  metal 
production  soon  to  follow.  Three  of  the  world's  largest 
brass  firms  were  located  on  its  banks  and,  at  the  turn  of 
the  century,  sixty  percent  of  the  world's  brass  products 
produced  annually  were  manufactured  in  the  valley.  As  a 
result  of  the  light  metal  production  industry,  the 
electroplating  industry  emerged,  and  approximately  30,300 
cubic  meters/day  of  treated  electroplating  waste  is 
currently  being  discharged  into  the  Naugatuck  River. 

Most  industrial  discharges  on  the  river  currently 
adhere  to  Best  Available  Technology  (BAT)  guidelines. 
However,  both  long-term  biological  monitoring  and  recent 
toxicity  testing  have  shown  toxic  impacts  of  metals.  As  a 
result,  the  Connecticut  Department  of  Environmental 
Protection  is  presently  developing  allocation  guidelines 
for  the  discharge  of  metals  into  the  river.  These 
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guidelines  will  be  based  upon  results  of  mathematical 
modeling  being  performed  by  the  State  Department  of 
Env i ronmenta 1  Protection.  The  data  set  available  was  used 
herein  to  calibrate  MEXAMS  to  the  Naugatuck  River  and  the 
output  will  provide  a  direct  comparison  of  two  model 
formulations  for  predicting  fate  and  transport  of  heavy 
metals  in  a  riverine  ecosystem. 


i 


l. 


CHAPTER  2 

MODEL  DESCRIPTIONS 
2.1  EXAMS  General  Description 

The  Exposure  Assessment  Model  System  (EXAMS)  was 
developed  in  response  to  a  need  to  describe  fate, 
transport  and  impacts  of  organic  contaminants  in  aquatic 
systems.  The  program  is  a  deterministic,  predictive 
systems  model,  based  on  a  core  of  mechanistic  process 
equations  derived  from  fundamental  theoretical  concepts 
(Burns  et  al . ,  1981 ) . 

The  model  estimates  exposure,  fate,  and  persis¬ 
tence  following  a  priority  pollutant  metal  discharge. 
Exposure  is  defined  as  the  steady-state  pollutant  con¬ 
centration  after  loadings  are  balanced  with  transport  and 
transformation  processes.  Fate  is  defined  as  the  dis¬ 
tribution  of  the  pollutant  in  the  system,  and  persistence 
is  the  time  required  for  effective  purification  of  the 
system.  The  program  uses  conservation  of  mass  to  balance 
loadings,  transport  and  transformation  of  the  compound. 

The  chemical  loadings  can  be  entered  into  EXAMS 
by  way  of  : 

l.j  Point  Source  Loadings 


—  ml*  y  i  n  j 
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2.)  Non-Point  Source  Loadings 

3.,  Aerial  Drift 

4. )  Atmospheric  Washout 

5. ;  Groundwater  Seepage 

These  chemicals  can  then  be  transformed  by  photolysis, 
hydrolysis,  biochemical  and  oxidation  reactions  in  par¬ 
ticular  compartments  of  the  environment  under  simulation. 
In  addition  to  the  above  transformation  processes,  EXAMS 
is  capable  of  three  transport  processes:  advection 
(dependent  on  the  hydrologic  data  input;,  dispersion 
(dependent  on  the  description  of  each  dispersive  pathway) 
and  volatilization  (the  air/water  interface;.  Advective 
and  dispersive  transport  are  possible  for  dissolved 
species,  sediment  sorbed  material  and  bio-sorbed 
materials.  Transport  through  the  ecosystem  compartments 
may  be  represented  by  whole  sediment  bed  loads,  suspended 
sediment  washloads,  exchanges  with  fixed  volume  sediment 
beds  and  ground  water  infiltration. 

EXAMS  can  also  be  used  to  simulate  the  ionization 
and  sorption  of  up  to  fifteen  molecular  species  of  an 
organic  chemical.  The  chemical  can  be  the  parent  un¬ 
charged,  singly  or  doubly  charged  cations  and  anions, 
each  of  which  can  occur  in  the  dissolved,  sediment-sorbed 
or  bio-sorbed  forms.  However,  EXAMS,  like  all  models, 
contains  important  assumptions.  These  assumptions  are 
listed  as  follows: 


l.>  A  first  order  evaluation  can  be  executed 


Figure 


independently  of  the  chemicals  actual  effects 
on  the  system, 

2. )  Linear  sorption  isotherms  and  second  order 

expressions  are  used,  rather  than  Michaelis- 
Menten-Monod  expressions,  for  biolysis, 

3. )  EXAMS  was  developed  for  long-term  or  steady- 

state  loadings  rather  than  dynamic  transitory 
(spills)  loadings,  and 

4. ;  Sorption  is  treated  as  an  equilibrium 

process,  that  is,  sorption/desorption  are 
assumed  to  be  rapid  compared  to  other 
processes. 


EXAMS  is  capable  of  simulating  both  lotic  and 
lentic  aquatic  systems.  These  systems  are  initially 
compartmentalized  and  classified  as  littoral,  epilimnion, 
hypolimnion  or  benthic  (Figure  1  illustrates  compartmen- 
talization  of  a  lake  system  and  a  river  system).  Each 
compartment  is  assumed  to  be  completely  mixed.  The  com¬ 
puter  then  develops  single  differential  equations  for 
each  compartment.  The  form  of  the  equation  and  the  vari¬ 
ables  are  defined  as: 


Vlf)  -  Le  +  Li  ’  V*C 


where: 

V  =  volume  of  water  in  compartment  (1) 

C  =  total  chemical  concentration  in  compartment  (mg/1) 
L  =  total  external  loading  (mg  hr; 

=  total  internal  loading  (mg  hr) 

K  =  overall  pseudo-first  order  loss  constant  (/hr) 
t  =  time 


EXAMS  was  not  originally  designed  to  describe 
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metal  transport  and  fate  and  has  been  modified  so  that  it 
can  interface  with  MINTEQ.  The  first  modification  in¬ 
volved  bypassing  calculations  that  are  unnecessary  for 
heavy  metal  speciation  or  duplicated  by  MINTEQ.  The 
second  modification  included  provisions  for  describing 
migration  and  fate  of  precipitated  fractions.  Ionization 
is  not  utilized  in  EXAMS,  because  MINTEQ  calculates 
dissolved  species  in  each  compartment.  Finally,  priority 
metals  are  generally  unaffected  by  the  transformation  and 
degradation  processes  that  effect  the  fate  of  organics  in 
natural  waters  and  these  processes  are  bypassed  in  the 
modified  EXAMS  version. 

2.2  MINTEQ  General  Description 

MINTEQ  is  a  computer  program  for  calculating 
aqueous  geochemical  equilibria  and  was  developed  for 
incorporation  into  MEXAMS.  The  major  objective  of  MINTEQ 
is  the  development  of  a  predictive  methodology  for  the 
assessment  of  the  metal  speciation  including  sorption, 
precipitation  and  ion  exchange  of  priority  pollutant 
metals  in  aquatic  environments.  This  chemical  equilibrium 
problem  is  described  as  a  set  of  mass  balance  equations, 
one  for  each  component,  and  a  set  of  mass  action  equa¬ 
tions,  one  for  each  specie.  The  equilibrium  constant 
approach  is  utilized  to  solve  the  equilibrium  problem. 
This  method  solves  nonlinear  mass  action  expressions 


using  linear  mass  balance  equations.  The  typical  mass 
action  and  mass  balance  equations  are  as  follows: 


Mass  Action: 


Vici  =  *i  I  X,  a(i, j) 
j=l  3 


(2) 


where: 

Y,  =  activity  coefficient  of  species  i 
c]  =  concentration  of  species  i 
1C  =  equilibrium  constant  reaction  i 
=  activity  of  component  j 
=  stoichiometric  coefficient  of 
components  in  species  i 
n  =  number  of  components 


Mass  Balance: 


m 

T ■  -  T  a(i, j)  C.  (3) 

3  j=l  x 


where: 


T.  =  total  analytical  concentration  of 
3  component  j 
m  =  number  of  aqueous  species 


Both  of  the  above  equations  are  definitions  for  systems 
with  no  solid  phases.  Solid  phases  are  dealt  with  using 
the  "transformation  of  basis"  method  (Felmy  et  al., 
1983).  This  method  reduces  the  number  of  independent 
variables  to  be  determined  and  allows  the  solution  of  a 
wider  range  of  chemical  equilibrium  problems. 


MINTEQ  is  also  capable  of  modeling  adsorption  in 
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a  number  of  ways.  Six  different  algorithms  are  accessible 
through  the  use  of  HI NTEQ  for  describing  sorption 
phenomena: 


ngmiur  equation 

3. )  "activity"  Freundlich  equation 

4. )  ion  exchange 

5. )  constant  capacitance  surface 

complexation  model 

6. )  triple  layer  surface 

complexation  model 


1 .  ;  "activity"  K  . 

2. )  "activity"  L§ 


The  "activity"  approach  is  the  ratio  of  the  amount  of 
metal  adsorbed  to  the  amount  remaining  in  solution  (ac¬ 
tivity  of  the  free  ion).  The  equation  can  then  be 
written  as  follows: 


SCu 

{cGy 


(4) 


where: 


K.  =  distribution  coefficient 
SCu  =  particulate  concentration  (molal  units) 
{ C u }  =  activity  concentration  of  the  free 
metal  (molal  units) 


The  "activity"  Langmuir  equation  has  an  advantage  over 
the  approach,  because  a  mass  balance  on  surface  sites 
is  considered.  The  "activity"  Langmuir  equation  is  writ¬ 
ten  as  follows: 


Kl  st  (Ol} 

SCU  =  — - - - 

1  +  {Cu} 


(5) 


Sj  =  maximum  quantity  that  can  be  adsorbed 
=  Langmuir  adsorption  constant 

The  third  approach,  the  "activity"  Freundlich,  is  used 
when  laboratory  data  obtained  does  not  conform  to  the 
"activity"  Langmuir.  The  "activity"  Freundlich  equation 
is  written  as  follows: 

1/n 

SCu  =  Kf  (6) 

where: 

K  r  -  Freundlich  adsorption  constant 
n  =  a  constant 

The  above  three  activity  approaches  are  preferred  over 
their  concentration  counterparts,  due  to  their  ability  to 
deal  with  a  wider  range  of  natural  waters.  However,  the 
activity  approach  is  not  without  limitations.  The  ac¬ 
tivity  approach  does  not  charge  balance  surface  sites  and 
adsorbed  species.  Secondly,  electrostatic  interactions 
between  adsorbing  ion  and  charged  surfaces  are  not  con¬ 
sidered.  Finally,  the  reaction  of  solid  with  aqueous 
constituents  other  then  adsorbate  ion  is  not  considered. 
In  order  to  reduce  the  above  limitations,  surface  com- 
plexation  models  have  been  developed  which  may  be  more 
applicable,  but  also  more  complex  (Felmy  et  al.,  1983). 
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The  constant  capacitance  model  and  triple  layer 
models,  both  based  on  theoretical  approaches,  consider 
charge  potential  differences  and  changing  properties  of 
the  surface  with  changing  pH  and  ionic  strength.  Both  of 
these  methods  require  a  great  deal  of  detailed  laboratory 
investigation  to  determine  the  critical  parameters  and 
are  primarily  research  tools. 

The  ion  exchange  approach  is  only  useful  when  the 
selectivity  coefficients  for  the  exchange  reactions  are 
known.  These  coefficients  are  generally  only  available 
for  bulk  electrolyte  ions  such  as  sodium,  potassium  and 
calcium.  The  ion  exchange  model  is  inadequate  for 
predicting  metal  adsorption,  because  metals  usually  form 
covalent  bonds  with  surface  sites  (Felmy  et  al.,  1983,. 

Finally,  the  equilibrium  composition  of  an 
aquatic  environment  is  determined  by  minimization  of 
Gibbs  free  energy  of  the  system  within  the  mass  balance 
constraints.  This  chemical  equilibrium  will  determine 
metal  speciation  in  addition  to  precipitation/dissolution 
and  sorption,  and  is  used  in  EXAMS  to  determine  the  metal 
fate  and  migration. 

2.3  MEXAMS  General  Description 

The  two  previous  sections  described  the  two 
individual  models  that  have  been  linked  to  form  the 
MEXAMS  model  (Figure  2).  The  MEXAMS  program  can  be  used 


mexams  flowchart 


Figure 
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in  both  batch  and  interactive  mode  with  three  different 
modes  of  operation  available.  The  first  mode  is  the 
MINTEQ  only,  which  allows  the  operator  to  determine  how 
changes  in  water  chemistry  will  affect  metal  speciation 
and  solid  phase  interactions  without  regard  to  transport 
processes.  Secondly,  the  EXAMS  only  mode,  which  deals 
with  the  ionization,  sorption,  transport  and  transforma¬ 
tion  of  a  given  pollutant.  Finally,  the  third  mode,  which 
links  EXAMS  and  MINTEQ  together,  allows  the  user  to 
determine  the  effect  of  transport  processes  and  chemical 
interactions  on  priority  pollutant  concentration.  This 
section  will  describe  the  operation  and  use  of  the 
coupled  mode. 

The  first  step  is  to  create  an  EXAMS  input  file 
that  describes  the  characteristics  of  the  aquatic  en¬ 
vironment  being  assessed.  The  user  then  enters  the  MEXAMS 
Interactive  Software  Package  (MISP).  MISP  will  call  for  a 
MINTEQ  file  for  each  compartment  that  contains  different 
water  quality  data.  The  programmer  will  also  input  run- 
specific  information  that  controls  the  number  of  times 
MINTEQ  updates  metal  concentrations.  MEXAMS  is  now  set  to 
simulate  metal  behavior,  migration  and  fate. 

Initial  calculations  by  EXAMS  are  made  assuming 
no  adsorption  or  precipitation,  and  the  initial  distribu¬ 
tion  of  dissolved  metals  in  each  compartment  is  obtained. 
These  concentrations,  along  with  suspended  sediment  and 


biota  concentrations,  are  sent  on  to  MISP.  MISP  then 
passes  the  dissolved  metal  concentration  in  the  first 
compartment  on  to  MINTEQ,  which  simultaneously  reads  the 
water  quality  data  for  that  compartment.  MINTEQ  then 
improves  the  estimate  of  the  concentration  of  metal  in 
solution,  adsorbed  to  sediments  and  biota  and  in  a 
precipitated  form.  These  results  are  returned  to  MISP  and 
metal  fractions  are  determined. 

MISP  then  proceeds  to  the  second  compartment  and, 
unless  the  water  quality  is  different  and  the  total  metal 
varies  by  more  than  5»,  MINTEQ  is  not  called.  MISP  util¬ 
izes  the  same  fractions  from  the  previous  compartment.  If 
total  metal  varies  by  more  than  5%,  MINTEQ  is  called,  but 
an  entire  new  calculation  is  only  performed  when  the 
water  quality  data  differs  from  one  compartment  to  the 
next. 

After  metal  fractions  have  been  computed  for  all 
compartments,  MISP  sends  the  information  to  EXAMS  and  the 
metal  transport  and  transformation  are  computed.  EXAMS 
then  predicts  new  exposure  levels  and  iterates  with 
MINTEQ  again.  Calculations  of  fate  and  persistence  are 
not  performed  until  all  iterations  are  complete. 

EXAMS  then  sums  the  flux  rate  of  metal  at¬ 
tributable  to  each  transport  and  transformation  process, 
in  which  the  significance  of  each  process  is  then  ob¬ 
tained  by  dividing  flux  rate  by  the  total  of  the  external 
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loadings.  Finally,  the  persistence  is  calculated  by 
terminating  the  external  loading  and  computing  the  dis¬ 
sipation  over  two  system  half  lives.  In  order  to  reduce 
computer  time  the  two  half  lives  have  been  divided  into 
twelve  equal  time  increments.  Since  metal  concentration 
will  change  with  time,  the  updating  of  each  compartment 
metal  fraction  is  possible.  The  user  specifies  the  number 
of  times  the  metal  fractions  are  to  be  updated.  The  metal 
fractions  will  change  with  each  time  increment,  and  thus, 
after  each  persistence  calculation  they  should  be  up¬ 
dated.  The  EXAMS  output  provides  exposure,  fate  and 
persistence  estimations,  the  MINTEQ  output  provides 
chemical  interaction  information  and  the  MISP  output 
provides  information  on  EXAMS-MINTEQ  interactions. 

The  MEXAMS  computer  program  is  capable  of  solving 
a  wide  variety  of  pollution  problems.  A  list  of  possible 
uses  follows: 


1. >  Rapid  evaluation  of  synthetic  organic 

pol lutants 

2. )  Rapid  hazard  evaluation  for  priority 

pollutant  metals 

3. j  Assess  impact  of  point  source  discharge 

from  mine  drainage 

4. ;  Provide  support  interpretations  of  metal 

bioassays 

5. }  Framework  for  what  is  or  what  is  not  known 

about  the  behavior  of  priority  pollutant 
metals  in  aquatic  systems 

6.  )  Management  tool 


Now  that  the  three  individual  programs  have  been 


presented,  the  specific  input  files  used  for  the 
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Naugatuck  River  will  be  discussed.  This  discussion  will 
include  a  description  of  the  compartmentalized  system 
used,  and  a  line  by  line  explanation  of  the  computer 
inputs.  The  explanation  will  state  any  limitations,  the 
variable  names  and  their  computer  format  code.  Some  of 
the  parameters  were  not  used  in  this  report  due  to  the 
fact  that  the  pollutants  are  priority  metals,  as  opposed 
to  synthetic  organic  chemical  contaminants. 


CHAPTER  3 


MODEL  INPUTS  AND  OUTPUTS 
3.1  Input  Description  (MINTEQj 

This  section  will  present  information  regarding 
the  operation  of  the  MINTED  program.  The  MINTEQ  program 
has  a  great  deal  of  flexibility  in  the  way  the  user 
defines  the  chemistry  of  the  ecosystem  being  modeled.  The 
program  is  capable  of  dealing  with  six  different  specie 
types,  as  described  in  the  MINTEQ  manual.  The  Type  I 
specie  is  the  component.  The  component  is  defined  as  the 
chemical  specie  chosen  to  represent  each  chemical  con¬ 
stituent.  Each  component  corresponds  to  a  MINTEQ  iden¬ 
tification  (ID)  number.  Table  1  provides  a  complete  list 
of  the  components  in  MINTEQ.  Type  II  complexes  are 
defined  as  aqueous  species  which  are  combinations  of  two 
or  more  components.  These  complexes  are  not  input  by  the 
user,  but  instead  are  output  by  the  computer  as  a  result 
of  the  Type  I  species  input.  Type  III  species  are  defined 
as  a  species  with  a  fixed  activity. 

Some  examples  of  Type  III  species  are  components 
at  fixed  activity  such  as  pH  or  pE,  solid  phases  present 
in  infinite  supply,  gases  present  at  a  fixed  partial 


pressure,  and  redox  reactions  between  two  components. 

Type  IV  species  are  precipitated  solids  subject  to  com¬ 
plete  dissolution.  The  Type  III  and  Type  IV  solids  are 
similar,  except  that  Type  IV  solids  have  a  finite  mass. 
Type  V  species  are  dissolved  solids  subject  to  precipita¬ 
tion.  These  solid  phases  will  only  precipitate  when  the 
saturation  index  is  exceeded.  Finally,  Type  VI  species 
are  species  not  considered  during  equilibrium  calcula¬ 
tions.  Instead,  these  species  are  considered  after  the 
equilibrium  calculations  are  complete. 

An  example  input  file  for  the  Naugatuck  River  is 
presented  in  Table  2,  with  the  numbers  on  the  left  cor¬ 
responding  to  line  numbers,  for  explanation  purposes.  The 
first  two  lines  are  descriptive  information  for  the 
sample  site.  The  third  line  describes  the  water  tempera¬ 
ture  iCj,  analytical  units,  and  fixed  ionic  strength 
designation.  If  ionic  strength  is  not  to  be  fixed,  0.00 
is  input.  Line  four  contains  nine  run  specific  user 
options,  inorganic  carbon  input,  debug  print  option, 
charge  balance  option,  considered  solids  and  print  op¬ 
tion,  number  of  iterations,  pH  variation  option,  fixed 
ionic  strength,  numerical  method  selection  and  output 
option.  The  fifth  line  contains  adsorption  paramaters. 
Line  six  through  24  are  the  component  input  lines,  and 
the  three  fields  correspond  to  the  specie  ID  number, 
total  analytical  mass  and  a  guess  for  the  log  of  the 
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TABLE  2 


MINTEQ  INPUT 


1.  NAUGATUCK  RIVER  CONDITIONS  (16-17  AUGUST  1983) 
COMPARTMENT  1  2N-CU-PB  ADSORPTION 

2.  KD  MODEL:  MATRIX  OBTAINED  FROM  BEST  AVAILABLE 
DATA  PH=  7.14 

3.  25.00  MG/L  0.00 


4. 

0  0  10 

0  0  0  0 

5. 

1  o.oc 

1  0.00 

0.00 

6. 

150 

2 . 000E+01 

-3.80 

7. 

460 

3 . 500E  +  00 

-3.90 

8. 

500 

1 . 890E+01 

-2.50 

9. 

140 

3 . 422E  +  01 

-5.50 

10. 

732 

2 . 500E  +  01 

-2.90 

11. 

180 

3.800E+01 

-2.80 

12. 

090 

8. 361E-02 

-5.10 

13. 

410 

3 . 090E+00 

-4.25 

14. 

580 

4 . 906E-01 

-11.0 

15. 

770 

8. 320E+00 

-6.50 

16. 

231 

5 . 450E-02 

-9.00 

17. 

950 

1 . 750E-02 

-8.00 

18. 

492 

7. 083E+Q0 

-8.00 

19. 

470 

3 . 550E-02 

-7.00 

20. 

281 

2 . 140E-01 

-7.00 

21. 

600 

6 . 0QQE-03 

-8.00 

22. 

330 

O.OOOE-OO 

-7.14 

23. 

991 

0 .  OOQE-OO 

0.00 

24. 

990 

0 .  OOOE  +  OO 

0.00 

25. 

26. 

3  3 

27. 

330 

7.14 

28. 

990 

29. 

991 

30. 

31.  2  3 


32. 

33. 

9909500 

2 

SOH-ZN 

1.00 

950 

1.00 

990 

-.24 

34. 

35. 

36. 

9916000 

2 

SOH-PB 

1.00 

600 

1.00 

991 

1.33 

37. 

38. 

39. 

9912310 

2 

SOH-CU 

1.00 

231 

1.00 

991 

0.20 

40. 

41. 
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component  activity.  This  section  is  the  water  quality 
data  for  each  compartment  of  the  ecosystem.  The  water 
quality  data  for  the  Naugatuck  River  was  obtained  from 
three  sources  ^A.J.  Medine  et  al.,  1982,  Curtis  1982,  and 
USGS  Water  Resources  Oata-Connecti cut ,  1982).  The  water 
quality  data  used  is  presented  in  Table  3.  The  next  line 
is  a  blank  line  separating  specie  Type  I  from  specie  Type 
III.  Line  26  corresponds  to  the  specie  type  and  the 
number  of  species  of  that  type.  The  next  three  lines  are 
for  fixed  pH  and  two  surface  sites  for  adsorption.  After 
another  blank  line,  the  adsorption  information  is  input. 
This  input  section  consists  of  nine  lines.  The  first  line 
contains  the  type  of  modification  and  number  of  these 
modifications.  Then  each  adsorption  input  has  three 
corresponding  lines  of  input.  The  first  line  is  the 
adsorption  description  and  the  log  value  for  the 
reaction. 

The  adsorption  routine  chosen  in  this  case  was 
the  "activity"  approach.  The  log  values  were  ob¬ 

tained  by  running  the  MINTEQ  program  without  adsorption. 
By  using  the  activity  results  and  the  field-observed 
particulate  concentration  for  each  metal  in  Equation  4, 
the  "activity"  was  obtained  ^log  was  then  utilized 

as  an  input  to  the  program).  Table  4  presents  the  data 
obtained  for  each  compartment,  and  the  respective 
values.  The  second  line  corresponds  to  the  components  and 
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Table  3 


Water  Qua!  ity  Data 


Concentrat ion 

Component  { mg / 1 )  Data  Source 


150 

- 

Ca 

1.50  E  +  01 

460 

- 

Mg 

3.50  E  +  00 

500 

- 

Na 

1.89  E+01 

410 

- 

K 

3.09  E+00 

470 

- 

Mn 

* 

281 

- 

Fe 

* 

770 

- 

S  i  02 

8.32  E+00 

950 

- 

Zn 

* 

090 

- 

B 

8.36  E-02 

231 

- 

Cu 

* 

732 

- 

S  i  0  4 

2.50  E+01 

180 

- 

Cl 

3.80  E+01 

140 

- 

C03 

3.42  E+01 

492 

- 

N03 

7.08  E+00 

580 

P04 

4.91  E-01 

600 

Pb 

* 

330 

- 

H 

* 

U.S.G.S.  (1983) 
U.S.G.S.  (1983) 

Medine  &  Conway  (1982) 
Medine  &  Conway  v 1982) 
Conn.  -  DEP  (1982)* ** 
Conn.  -  DEP  (1982)** 
U.S.G.S.  (1983) 

Conn.  -  DEP  (1982;** 
Medine  &  Conway  (1982) 
Conn.  -  DEP  (1982)** 
U.S.G.S.  (1982) 
U.S.G.S.  (1982; 

Medine  4  Conway  (1982) 
U.S.G.S.  (1982; 
U.S.G.S.  (1982; 

Conn.  -  DEP** 

Conn.  -  DEP** 


*These  inputs  are  compartment  dependent 

**Connecticut  State  Department  of  Environmental 
Protection 


» 
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TABLE 

4 

* 

1 

Compartment 

Kd  Data 

r 

Compart- 

U 

ment 

Free  Ion 

Particulate 

#-Type 

Element 

Activity  (H) 

Cone.  (M) 

Log  Kd 

1-E 

Cu 

9.741  E - 0  8 

1.50  E-07 

.20 

► 

Pb 

9.940  E-10 

2.17  E-08 

1.33 

Zn 

1.321  E - 0  7 

7.65  E-08 

-.24 

. 

2-B 

Cu 

8.467  E - 0 7 

4.41  E-03 

3.73 

:• 

Pb 

9.300  E - 0 7 

2.79  E-03 

3.48 

Zn 

5.451  E - 0  5 

3.89  E-03 

1.85 

B 

3-E 

Cu 

1.634  £-07 

1.34  E-07 

-.09 

- - :  — 

Pb 

2.551  E - 0  9 

2.17  E-08 

.93 

Zn 

1.830  E-07 

9.94  E-08 

-.26 

4-B 

Cu 

1.796  E-06 

1.59  E-03 

2.95 

to 

Pb 

Zn 

1.834  E-06 
4.637  E-05 

3.00  E - 0 4 
3.49  E-03 

2.21 

1.88 

«*  '-‘--•to 

5-E 

Cu 

4.422  E - 0 8 

3.93  E-08 

-.05 

Pb 

2.234  E -09 

2.17  E-08 

.99 

Zn 

9.065  E - 0 8 

5.35  E-08 

-.23 

• 

6-8 

Cu 

3.225  E-06 

1.51  E-03 

2.67 

Pb 

2.543  E-06 

3.07  E - 0 4 

2.08 

* 

Zn 

2.177  E-05 

3.66  E-03 

2.23 

• 

• 

E  -  E pi  imn ion 

B  -  Benthic 

i 

i 

the  stoichiometry  of  the  reaction.  Finally,  the  next  two 
lines  are  blank  to  signify  the  end  of  the  input  file. 


This  section  described  the  input  file  for  the 
Naugatuck  River.  Also  included  in  the  appendices  is  a 
copy  of  the  MINTEQ  input  description  (Felmy  et  al., 

1983; .  This  description  includes  computer  format,  pos¬ 
sible  options  and  a  more  in-depth  description  of  the 
MINTED  input  file  itself.  Following  tne  input  of  tne 
MINTEQ  data  file,  the  next  step  is  the  creation  of  an 
EXAMS  input  file. 

3.2  Input  Description  ^EXAMSj 

The  EXAMS  input  file  is  generated  based  on  a 
compartmental ization  of  the  riverine  ecosystem.  EXAMS,  as 
stated  previously,  is  responsible  for  the  metal  fate  and 
migration  calculations.  Figure  3  displays  the  initial 
compartmental ization  scheme  used  for  the  Naugatuck  River 
system.  Once  the  compartmental ization  is  selected,  the 
EXAMS  input  file  is  generated.  The  file  consists  of  67 
lines  of  input,  which  will  be  explained  next  (Table  5). 

The  first  two  lines  are  descriptive  lines  for  the 
ecosystem,  the  chemical  name  and  ecosystem  name,  both  a 
maximum  of  three  characters.  The  third  line  corresponds 
to  the  number  of  chemical  loadings.  In  this  case,  there 
are  five  loadings  and  the  next  five  lines  describe  the 
loading  as  one  of  five  types  (units  of  input  for  each 
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Figure 
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Table  5 


EXAMS  Input 


MAG 


1 . 2646 -02 
0. 

4.4046-03 
3.  361E-G2 
C. 


0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 


0. 

0. 

0. 

0. 

0. 


6. 505E-04 
0. 
0. 

3. 5386-04 


hvm 

He i V  y  Me t  6 ' 

C  0  0  0  c 

NAu 

ftAjGATUC*  RIvES 
6 
L 
8 
L 
8 
L 
6 

41  80 

G.OOGC  0.0000 

8  3160-02  4.9256- 

6  3666*04  5. 3860 ■ 
0.1706?  C.00814 

9  6666-020. 00 C 
2  9666  c .  or-: 

0  000  o.coo 

o.ocr  o.ooo 
1  6266*02  1.6266 


CALIBRATION  TEST  RUN 


0  0000  0. 0000  0. 0000  0. 0000 
02  3  3536-03  1.1496-02  7.9036-04  2.4026-03 
04  1  257E-04  1.2576-04  2.627E-05  2.6276*05 


7963  0.00914 

4  6946-01  0.000 

ooc 


0.  304e 


0.489 

0.  OOC 
0.  I/O 
2  C  OOC 


OOC 

o.  c :  g 
0.  ooc 


02.  00 
IOC.  0 
0.000 
0.000 
0.  000 
c.  ooc 
0,000 
oc  cc 
0 

:■  oc; 
o.ooo 

e 


c.ccc 
c.  coo 
0.  ooc 

«  r  rtf 
V  .  N.  V.  - 

0.  OC  0 

or. 
OC  00 
1.  532- 

115.2 


000 

002 

00? 

000 

oc. 


0.  coc 
0.  000 
c.  ooc 
0.  ooc 
0.  000 


oc 


3  OCO 

o.ooc 


oc  cc 
00.  02 


1 .  530 
lii 


00.  00 
00.  oc 


0.00914 
7.  7086*02  O.OOC 
5.528  0.000 

0.000  0.000 

O.OOC  0.000 

3.5386-02  3.5386-02 


0.  000 
0.000 
0.000 
0.000 
o.ooc 


0.  000 

0.  ooc 

o.  occ 

0.000 

0.000 

0.000 

0.000 

0.  000 

0.000 

o.ooc 

0.000 

0.000 

0.000 

0.000 

00.  00 
00.00 


7.  SCO 
0.000 


.  530 
115.2 


00.00 

00.00 


1  3 

3  5 

1  000 

3 

5  2 

0  1 

1 .  ooc 

4  4 

3  6 

1  000 

6 

5 

0.900 

2 

4 

0.  900 

0.100 

J  3 

2  4 

i  Sbbi+OA 
b.  992 1 - 02 
i.  000 £-07 
36C0E-01C 

5 

6 

1.2576*04 

5.4496-02 

5.0006-07 

'.38006-010 

2. 6276-04 

1 . 5706-01 
5. 000E-C7 
. 38006  -  01 C 

. 38006- 

■010. 38006 

-010. 3800E-01 

2b.  00 

25.  00 

25.00 

25.00 

25.  00 

25.  00 

2b  00 

25.  00 

25.  00 

25.  00 

25.  00 

25.  00 

2S.  SO 

23.  50 

2  3.  50 

23.50 

23.  50 

23.50 

7.143 

7.  142 

7.007 

7.007 

7.078 

7.078 

e.  8S7 

6.85' 

6.993 

6.993 

6. 922 

6.922 

3. 1 000  £ -08 

0. 00. 10006-080. 5000E 

-09 

0.00. 10006-08 

OC.  00 

00.00 

00.00 

00.00 

00.00 

00. 00 

OC  .  00 

00.00 

00.00 

00.00 

00.  00 

00.  OD 

23 .  SO 

23.50 

23.50 

23.50 

23.50 

23.50 

00 .  OC 

00.  00 

00.00 

00.00 

00.00 

00.  00 

OC .  CO 

00.00 

00.00 

00.00 

00.00 

00.  00 

00 .  00 

GO.  00 

00.  00 

00.00 

00.  00 

00.  00 

00.00 

00.  OC 

00.00 

00.00 

00.00 

00.  00 

00 . 00 
00.00 

00.  OC 

00.00 

00.00 

00.00 

00. 00 

00.  00 

00.  00 

00.00 

00.00 

00. 00 

1 .  ooc 


0.102 


28 


chemical  loading  are  Kg/hry.  The  five  types  of  loadings 
are  point  source,  non-point  source,  rainfall,  aerial 
drift  and  ground  water  seepage.  The  loadings  are  input 
from  left  to  right  for  each  compartment  in  the  order 
previously  stated.  The  next  two  lines  are  a  description 
of  the  chemical,  followed  by  a  line  of  flags  denoting 
which  species  of  the  organic  toxicant  can  occur.  A  zero 
is  used  for  no  occurrence  and  a  one  for  occurrence,  and 
in  this  case,  zeros  are  input  because  the  pollutant  being 
investigated  is  a  heavy  metal.  The  next  two  lines  are  an 
ecosystem  description. 

The  ecosystem  was  divided  into  unequal  compart¬ 
ment  lengths  representing  well-defined  segments  based  on 
actual  point  loadings.  Compartments  one  and  two  are  1.3 
miles  long,  compartments  three  and  four  are  .3  miles  long 
and  compartments  five  and  six  are  6.7  miles  long.  Figure 
4  displays  the  four  different  sample  sites  (East  Albert 
St.,  Gulf  Stream,  Bogue  Road,  Bridge  Abutment)  within  the 
stream  segment  selected  for  compartmental ization  and 
model  i ng . 

The  next  seven  lines  are  the  number  of  compart¬ 
ments,  and  the  compartment  code  for  each  compartment.  The 
compartments  are  classified  as  epilimnion  (Ej,  hypolim- 
nion  (H),  littoral  (l)  or  benthic  (B).  The  compartment 
classification  is  then  followed  by  the  latitude  input  in 
degrees,  and  the  wind  velocity  at  10  cm  above  the  water 
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surface  tm/sec,  for  each  compartment.  Line  24  is  the 
environmental  volume  (cubic  meters;  for  each  compartment 
followed  by  entries  of  area  (square  meters^  and  average 
depth  (m)  for  each  compartment. 

The  next  section,  lines  27  through  31,  represent 
flows  entering  the  ecosystem  compartments.  These  flows 
are  of  the  following  type:  stream  flow  ^cubic  meters/hrj, 
stream-born  sediment  load  (Kg/hr),  non-point  source  water 
flow  vcubic  meters/hrj,  non-point  source  sediment  loads 
(Kg/hr),  and  ground  water  seepage  flow  (cubic  meters/hr). 
The  next  input  is  average  rainfall  (mm/moj  in  the 
geographic  area,  then  the  average  cloudiness  in  tenths  of 
full  sky  cover  is  input.  From  this  point  on  in  the  data 
set,  many  parameters  are  input  as  zero  because  they  deal 
with  transformation  of  organic  toxicants  and  not  priority 
pollutant  heavy  metals. 

The  next  six  lines  are  inputs  for  average 
spectral  irradiance  immediately  below  the  water  surface. 
These  inputs  are  not  required  for  transf ormat i on  of  heavy 
metals.  Line  40  corresponds  to  the  evaporative  water 
losses  (mm/moj  from  each  compartment.  The  next  line  is 
where  suspended  sediment  concentrations  (mg/1)  are  input 
if  dealing  with  water  column  compartments  and  bulk  den¬ 
sity  (g/cc)  for  benthic  compartments.  These  inputs  are 
followed  by  percent  water  in  bottom  sediments  of  benthic 
compartments.  Water  column  compartments  have  inputs  of 


zero  for  this  variable. 

The  transport  mechanisms  of  advection  and  disper¬ 
sion  are  input  in  lines  43  -  52.  Advection  is  input 
first,  with  the  first  line  correspond  ing  to  the  number  of 
advective  pathways.  The  next  three  lines  are  the  source 
compartment,  the  receiving  compartment  and  the  proportion 
of  the  advective  flow.  The  dispersive  pathways  are  then 
input.  The  first  line  equates  to  the  number  of  dispersive 
pathways.  The  source  compartments  are  then  entered  fol¬ 
lowed  by  the  receiving  compartment  and  the  description  of 
the  exchange  pathway.  The  cross-sectional  area  (square 
meters],  c haracter i st i c  length  and  eddy  dispersion 
coefficients  (square  meters/hr)  are  the  parameters  re¬ 
quired  to  describe  the  dispersive  pathway.  The  charac¬ 
teristic  length  is  defined  as  the  distance  from  compart¬ 
ment  center  to  compartment  center. 

The  next  three  lines  correspond  to  the  organic 
carbon  content  of  compartment  sediment,  cation  and  anion 
exchange  capacity  of  the  sediments.  Line  56  is  the 
average  temperature  for  each  compartment.  The  pH  and  pOH 
of  each  compartment  are  the  next  two  inputs.  Finally,  the 
remaining  lines  are  not  required  for  the  modeling  of  the 
Naugatuck  riverine  ecosystem  for  the  migration  and  fate 
of  priority  pollutant  metals. 

Included  in  the  appendices  is  a  copy  of  the  EXAMS 
input  file  description  (Burns  et  al.,  1981).  The  descrip- 
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tion  provides  a  detailed  description  of  the  computer 
format,  and  the  parameters  which  were  not  used,  but  can 
be  used  for  organic  toxicants. 

For  each  compartment  there  is  a  corresponding 
MINTEQ  input  file,  and  it  should  be  re-emphasized  that 
each  compartment  is  assumed  to  be  completely  mixed  after 
the  priority  pollutant  enters  the  compartment.  After 
input  files  have  been  created  for  each  compartment 
(MINTEQ),  and  the  transformation  and  transport  model 
(EXAMS)  input  file  is  complete,  the  MISP  interactive 
input  file  is  created. 

This  file  encompasses  twenty-five  lines  in  this 
example.  These  lines  allow  the  interaction  of  MINTEQ  with 
EXAMS  and  produce  the  MEXAMS  system  model.  The  first  line 
corresponds  to  the  selection  of  the  computer  mode  (MINTEQ 
only,  or  coupled  version).  The  metal  selection  is  made 
next  from  the  list  of  metals  (As,  Cd,  Cu,  Pb,  Ni,  Ag  and 
Zn)  and  corresponding  identification  numbers  in  the 
MINTEQ  components  list.  The  next  line  corresponds  to  how 
many  MINTEQ  files  are  needed.  The  number  of  required 
files  is  dependent  on  the  water  quality  data  for  the 
aquatic  system.  If  the  water  quality  data  does  not  change 
significantly  for  each  compartment,  then  one  file  is 
sufficient.  Water  quality  fluctuations  could  be  a  result 
of  wastewater  treatment  plant  discharges,  industrial 
discharges,  overland  runoff,  etc. 


The  next  group  of  lines  corresponds  to  the  number 
of  compartments  that  match  up  with  each  MINTEQ  file, 
followed  by  the  compartment  number  that  corresponds  to 
each  MINTEQ  file.  The  next  input  is  the  number  of  times 
the  programmer  wants  the  steady  state  metal  concentration 
updated.  EXAMS  initially  assumes  all  metal  is  in  the 
dissolved  form  and  by  providing  MINTEQ  with  these 
results,  the  actual  dissolved,  sediment  sorbed,  biosorbed 
and  precipitated  metal  concentrations  are  obtained.  EXAMS 
uses  these  updated  calculations  in  the  next  simulation 
for  fate  and  transport  of  the  pollutant.  Line  seventeen 
relates  to  how  often  MINTEQ  is  used  for  persistence 
calculations.  After  steady  state  has  been  achieved,  the 
pollutant  loading  is  ceased  and  persistence  calculations 
begin.  These  calculations  are  based  on  two  system  half- 
lives,  which  are  divided  into  twelve  equal  time  incre¬ 
ments.  This  input  will  continue  metal  fraction  updating 
throughout  the  persistence  calculations.  The  next  three 
inputs  correspond  to  the  output  option,  the  MINTEQ  file 
names  for  each  compartment,  and  the  EXAMS  input  file 
name.  There  are  three  output  options  available:  the  full 
MINTEQ  output,  aqueous  species  distribution  and  all  mass 
totals,  and  all  mass  totals.  This  concludes  the  MISP 
interactive  input  file.  Table  6  is  a  sample  of  the  MISP 


Table  6 


MISP  Input 


3 

231 

6 

1 

1 

1 

2 

1 

3 
1 

4 
1 

5 
1 

6 

3 

4 
3 

NAUG1.DAT 
NAUG2.DAT 
NAUG3.DAT 
NAUG4.DAT 
NAUG5. DAT 
NAUG6.DAT 
EXAMS21.DAT 


The  MINTEQ  output  data  received  when  the  un¬ 
coupled  mode  is  operated  can  be  divided  into  seven  sec¬ 
tions.  The  first  section  is  a  duplication  of  the  first 
five  lines  of  the  input  file.  Section  two  corresponds  to 
the  initial  charge  balance  of  the  solution  being  modeled, 
prior  to  aqueous  speciation.  A  large  percent  difference 
can  be  indicative  that  one  or  more  major  cations  or 
anions  may  have  been  neglected  when  preparing  the  water 
quality  input  file.  However,  the  water  analysis  may  be 
entirely  correct  with  the  final  speciation  indicating  an 
acceptable  charge  balance.  Therefore,  initial  charge 
imbalance  should  be  evaluated  on  a  case  by  case  basis 
(section  2  is  presented  in  Table  7j. 

The  third  section  provides  an  iterative  output 
displaying  the  convergence  of  the  first  component  in  the 
input  file  to  its  actual  speciated  activity  (Table  7). 
Also  included  in  this  section  is  the  calculated  molality, 
activity.  Log  activity,  activity  coefficient,  and  new  Log 
K  for  all  input  components.  The  new  Log  K  is  the  effec¬ 
tive  equilibrium  constant  after  temperature  and  ionic 
strength  corrections.  This  section  can  be  used  to  enable 
the  user  to  make  a  better  guess  for  the  activity  in  the 
input  file. 

The  fourth  section  provides  the  output  data  for 
the  six  different  specie  types.  The  output  for  Type  I 
components  and  Type  II  complexes  are  similar  with  the 


Table 


'■V  n  '! 


36 


r 

b 


10 

Z 

C 


< 

o 


co 

z 

o 

z 

< 

o  \ 
o  - 
»  <o 
UJ  z 
o  o 
(0  - 
®  z 
•  < 
CVJ 


<0 
(0  z 

—  H 

z  < 

<  o 


£  * 
5  o 
co  ■ 
o 

«  - 

o  i 

O  UJ 

uj  n 

b-  O  n 

<  CVJ 

—  .  UJ 

CPJU 
UJ  z 
CL  n  UJ 

o  co  cr 
z  z  UJ 

D  O  U 

-  u 

UJ  <  Q 

o  o 

Z  K 

<  U  Z 

J  O  UJ 

<  o 

mzx 

5  uj 

UJ  CO  0. 


g 

5 


K  o  o>  o  o  «  « 

>  O  Ifi  N  •-  CM  lf> 
hO'J^OOft 
uoeoMOM- 
<  CO  O  O  rf  IL  ft 


8 


o  «  co  o  o  co 


^  io  n  ^  *  ® 

z  o  o  o  o  o  o 

X  »  •  «  •  »  » 

U  UJ  UJ  UJ  UJ  UJ  UJ 

o  ®  o  ®  n  cti 

u.  —  o  <o  in  co  cm 

L.DN^fflO>  V 

o  cvj  6  -  cvj  -  co 


-J  o  o  o  o  o  o 

Q  i  •  i  i  i  • 

tuiujuiujuu 

JO'CIOIC'C  3> 
<  0)  O  0)  (J1  3» 

K 

>  k 

-j 

o 

y> 


K  UJ 
3  Z 
O  < 
Z 

(0 

z 

o 


-  o  o  o  o  o  o 

K 

< 

ujui  —  ruco^r® 


z 

x  n  r*  n  ® 

u  o  o  o  o 

i  i  i  i 

U  UJ  UJ  UJ  UJ 

u  o  ^  <r  to 

~  cvj  co  ^  o 

O  ®  O  (0  ® 
CO  -  0> 

till 

©  0>  ©  O  * 
OOOOM 
_)  O  O  CVJ  o 

~  -  o  - 

2  •  ■  ■  • 


<©©©«- 
r  ®  k  n  to 

z  ©  o  ©  © 

3  n  ^  © 

O  9)  0)  V  « 

N  N  0)  N 


o 

CO 

© 

« 

<0 

© 

to 

<? 

® 

10 

*— 

co 

•— 

10 

*- 

CVJ 

* 

© 

— 

to 

<o 

o 

n 

6 

> 

K  5T  51  ^  N 

-  o  o  o  o 

>1111 

—  UJ  UJ  UJ  UJ 

H  N  N  ^  O 

o-onoi 

<  ®  -  N  © 

n  -  rv  cvj 

«J  5*  ^  *  N 

“  O  O  O  O 


8 


UJ  UJ  UJ  UJ 

o  to  to  —  cvj 

j  -  ©  -  Cl 

«  ®  co  cnj 

o  •  •  • 

5J  —  ®  CO 


(0  J 


81? 


O  O  O  ' 


UJ  UJ  UJ  UJ 
M  J  -  O  CVJ  O 

<  ©  *  CVi  O 

co  z  ©  ^  cvj  n 

z  <  •  •  •  • 

o  ^  —  ®  m 


< 

z 

UJ 


-  Ul 

z 

< 

••  z 

< 

K 

§ 


<0 

<  0  <  o 
oczu 


D  O  O  O  O 

o.  r>  o  o  v 

D  - 
C 


.  t 
.  i 


37 


following  headings:  calculated  molality,  activity.  Log 
activity,  activity  coefficient  lGAMMA),  new  Log  K,  and 
the  enthalpy  of  the  reaction  (DH).  Type  III  fixed  solids 
are  of  four  basic  categories,  redox  reactions,  components 
at  fixed  activity,  solid  phases  and  gases  at  a  fixed 
partial  pressure.  The  output  for  these  type  species 
consists  of  calculated  molality.  Log  molality,  new  Log  K 
and  the  DH.  In  our  case  the  Type  III  species  fall  into 
the  category  of  fixed  activities,  thus  any  intuitive 
information  from  this  section  cannot  be  obtained. 

However,  if  solid  phases  were  being  dealt  with,  a  posi¬ 
tive  calculated  molality  would  equal  the  mass  of  solid 
precipitated.  On  the  other  hand,  a  negative  calculated 
molality  would  equal  the  dissolved  mass  or  mass  added  to 
the  solution.  No  Type  IV  or  Type  V  species  were  permitted 
from  the  list  of  solids  which  were  supersaturated  as 
indicated  by  MINTEQ.  Chemical  precipitation  has  not  been 
observed  in  the  Naugatuck  River  although  it  cannot  be 
precluded.  Type  V  calculated  molality  would  correspond  to 
saturation  index  values.  Type  VI  has  two  different  inter¬ 
pretations  for  the  calculated  molality.  If  solid  phase, 
it  corresponds  to  saturation  index.  If  aqueous  species, 
it  corresponds  to  activity  of  specie.  A  sample  of  section 
four  output  is  displayed  in  Table  8. 

Section  five  presents  the  percentage  distribution 
of  Type  I  components  and  Table  9  is  an  example.  The  sixth 
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Table  8 


MINTEQ  Output 


SPECIES  TYPE  I 
ID  NAME 

150  CA 

460  no 

500  NA 
140  C03 

732  504 

160  CL 
90  M3B03 

410  K 
380  P04 

770  H4S104 

23)  CU*2 
950  ZH 
492  NO 3 
470  MN-2 
20)  EE*3 
600  PB 

SPECIES  TYPE  I 

•  -> 

1303330  CAOH  ♦ 
150U00  CAHC03  ♦ 
ISO’JO)  CAC03  AO 
150^320  CAS04  AO 
1503600  CAHPC4  AO 
1305801  CAP04  - 
1503802  CAM2P04  ♦ 
3001 400  NAC03  - 
300’ 401  NAHC03  AQ 
600 J 600  PBCL  ♦ 
6001801  PBCl2  AO 
6001802  PBCL3  - 
600 ’ 803  PBCl4  -2 
600 ’J 00  PB l 003)2-2 
6003300  PBOH  - 
6003301  PB ( OH ) 2  AO 
6003302  PB I  OH 1 3  * 
6003303  PB20H  *3 
6004920  PBN03  ♦ 
600’’320  P3S04  AQ 

6003304  PB3 1  OH  J  4  *2 
600)401  PBC03  AQ 
6003335  PB < OH) 4  >2 
6007321  PB 1 S04 ) 2 - 2 
SPECIES  TYPE 

J  D  NAME 

2  H20 

991  S0H2 
990  SOH 1 
330  H 

SPECIES  TYPE 

l 0  NAME 

3301403  C02(GAS> 


-  COMPONENTS 
CAlC  mol 

4  615E-04 

1  393E04 
0  21 1 E-04 
3  792E-07 

2  408E-04 
1  072E-03 
1  341 E -06 

7  692E '05 
1  940E- 1 1 

8  643E-05 
1  223E-07 
1  67 1 E -  07 
1  )42 E-04 
6  220E-07 
1  245E- 14 
1 . 286E-09 


activity 

3  01 7E  *  04 

1  .  107E-04 
7 . 734E  * 04 

2  996E07 
1  091 E  *  04 
1 . 006E-03 
1  342E-06 
7 . 423E - 05 
1  . 1 37E  *11 
0. 649E-05 
9  6S6E-06 
1  323E-07 
1  074E-04 

4  925E  *  07 
7  545E-15 
1 . 010E-09 


LOO  ACTVTY 
*3.41030 
'3. 95593 
'3. 11159 
-6  52340 
-3.72333 
-2.99650 
-5  07224 
-4  12945 


I  I  -  complexes 

r>4i  r  ».rt"CTv 

1  409E-09  1  329E-09 

1 . 944E-06  1 . 634E-06 

1.624E-07  1  625E-07 

1  469E-05  1  470Z-05 

3  8226-07  3  825E-07 

1 . 325E-06  1  249E-O0 

2.2O4E-O0  2.O70E-O0 

4  555E-00  4.295E-09 

2. 017E-07  2  01 0E-O7 

4  307E '  1  1  4  055E-U 

6.475E-14  6.479E-14 

5  490E-17  5. J76E-17 

3  106E-2O  2.503E-20 

5  0C6E-I2  3.959E-12 

2  888E-10  2.7192-10 

1  45SE-12  1  460E-12 

2  450E-16  2.314E-16 

1  0582-17  6  152E-16 

1  705E-12  1.606E-12 

1  074E-10  1.074E-10 

6  282E*23  4.937E-23 

5  257E-09  5  26 1 E - 09 

9  334E-2I  7  335C-21 

1.357E-13  1  066E-13 

III  -  PIXEO  SOL I  OS 

CALC  MOL  LOO  MOL 
-9247E-06  -5  03 

-1.751E-07  -6  75 

-7  61 3E-O0  -7.11 

-6.J25E-04  -319 


-0.034 
-6  757 
-7.  110 
-3  192 


BAMMA 
. 792606 
.794370 
.941939 
. 709951 
.705340 
.940572 
1 . 000697 
. 940572 


NEW  LOOK 
.  101 
.  100 
.026 
102 
.  105 
.027 
-  .000 
027 


CALC  MOL 
2  273E-03 


LOO  MOL 
-2.643 


-10  94412 

.506236 

.232 

-4 . 06305 

1  000697 

- .  ooo 

-7  01305 

.791776 

.  101 

-6  07842 

. 791776 

101 

-3  96886 

940473 

027 

-6  30750 

.791776 

.  101 

-  14 .  12233 

606033 

218 

-6.99549 

.  705013 

.  105 

n***M«\ 

wn  •  • 

-8  07633 

. 943302 

-12  573 

-5  73669 

943302 

1  1  370 

-6.76917 

1  000697 

3  152 

-4  83263 

1  000697 

2  309 

-6  41 742 

1  000697 

15  005 

-7  90342 

942750 

6  405 

-7  08242 

.942750 

20  906 

-0.36707 

942750 

1  294 

-6.69507 

1  000697 

10  060 

-10  39199 

. 94 1 520 

1 . 626 

-13  16649 

1 . 000697 

1 . 600 

-16  20599 

94 1 520 

1  725 

-19  60149 

.705013 

1  .405 

-It  40245 

. 705013 

10  745 

-9  56552 

941520 

■7  684 

-11  03554 

1 . 000697 

-17  120 

-15.63557 

.941520 

-20  034 

-17  21101 

.501390 

-6  124 

-1  1  79437 

. 94 1 520 

1.196 

-9  96882 

1  000697 

2  750 

-22  30658 

705013 

-23  775 

-0.27897 

1  000697 

7  240 

-20  13460 

. 705013 

-39  594 

-12.97214 

.765013 

3  575 

NEW  LOOK 

OH 

000 

0 

0 

0 

0 

0 

7  140 

0 

RED 

NEW  LOOK 

DM 

18  160 

-.530 

14  535 

1  790 
4  030 
1.470 
-  230 
3.100 

-1  . 120 
0.91  1 
0 

4  300 

1 .060 

2  170 

3  530 


ooooooooooo 


section  presents  the  amount  of  each  Type  I  component  in 
solution  ^aqueous  mass),  sorbed  (sorbed  mass)  and 
precipitated  (precipitated  mass).  The  final  charge 
balance  is  also  presented,  and  if  the  final  percent 
difference  after  aqueous  speciation  is  less  than  five 
percent,  the  run  is  acceptable.  A  sample  of  section  six 
output  is  presented  in  Table  9.  The  heavy  metal  under 
investigation  in  this  case  study  is  copper,  and  from 
Table  9  it  can  be  seen  that  14.3  percent  of  the  copper  is 
in  the  free  ion  form,  16.9%  bound  to  surface  sites  with 
the  remainder  as  inorganic  complexes. 

The  final  section  provides  the  saturation  indices 
for  all  minerals  and  solids  in  the  data  bae.  This  sec¬ 
tion  is  used  to  determine  if  the  saturation  index  of  a 
particular  mineral  or  solid  has  been  exceeded.  A  positive 
Log  SI  means  the  solid  is  oversaturated .  This  also  tells 
you  that  this  solid  could  precipitate,  and  thus,  if  your 
environment  favors  this  solid,  the  solid  could  be  input 
as  a  Type  IV  solid  and  allowed  to  precipitate.  Finally,  a 
line  will  be  printed  "NORMAL  TERMINATION  OF  MINTEQ"  and 
the  output  is  complete.  The  final  section  is  presented  in 
Table  10.  An  example  of  the  entire  MINTEQ  output  is 
presented  in  the  appendices. 

3.4  Output  iEXAMS ) 


The  EXAMS  output  obtained  when  the  coupled  mode 
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Table  9 

MINTEQ  Output 


PERCENTAGE  DISTRIBUTION  OP  COMPONENTS 


CA 


96 

5 

PERCENT 

BOUND 

IN 

SPECIES 

ISO 

CA 

2 

9 

PERCENT 

BOUND 

IN 

SPECIES 

1507320 

CAS04  AQ 

CU*2 

14 

3 

PERCENT 

BOUND 

IN 

SPECIES 

231 

CU+2 

17 

9 

PERCENT 

BOUND 

IN 

SPECIES 

9912310 

SOM - CU 

id 

2 

PERCENT 

BOUND 

IN 

SPECIES 

2311400 

CUC03  AQ 

i 

7 

PERCENT 

BOUHO 

IN 

SPECIES 

231 3300 

CUOH  ♦ 

44 

9 

PERCENT 

BOUND 

1  N 

SPECIES 

2313301 

CU(0H)2  AO 

2 

6 

PERCENT 

BOUND 

IN 

SPECIES 

231 1402 

CUHC03  ♦ 

ZN 

62 

4 

PERCENT 

BOUND 

IN 

SPECIES 

950 

ZN 

2 

2 

PERCENT 

BOUND 

1  N 

SPECIES 

9507320 

ZNS04  AQ 

2 

9 

PERCENT 

BOUND 

1  N 

SPECIES 

9501400 

ZNHC03  ♦ 

3 

0 

PERCENT 

COUND 

IN 

SPECIES 

9501401 

ZNC03  AQ 

28 

4 

PERCENT 

BOUND 

IN 

SPECIES 

9909500 

SOH - ZN 

D* 

name 

AQUEOUS  MASS 

SORBED  MASS 

PRECIPITATEO  MASS 

150 

CA 

4  987E-04 

0 

0 

460 

NO 

1  430E-04 

0 

0 

500 

NA 

8  221 E- 04 

0 

0 

140 

C03 

5  702E-04 

0 

0 

732 

SG4 

2. 601 E- 04 

0 

0 

1  80 

CL 

1 . 072E-03 

0 

0 

90 

M3B03 

1  352E-06 

0 

0 

410 

K 

7  302E-05 

0 

0 

580 

P04 

5  1  SUE -06 

0 

0 

770 

M4SI04 

8  658E-05 

0 

0 

231 

CO-2 

7.042E-D7 

1  533E-07 

0 

950 

ZN 

1 . 015E-07 

7  613E-08 

0 

49  2 

NO  3 

1  142E-04 

0 

0 

470 

MN»2 

6  U5eE-07 

0 

0 

281 

FE*3 

3  832E-06 

0 

0 

600 

PB 

7  360E-09 

2  160E-08 

0 

2 

H20 

9  247E06 

0 

0 

330  M 

chapge  dalance 

6  425E -04 
SPEC) ATED 

0 

0 

Sun  OF  CATIONS  *  2  149E-03  SUJ1  OF  ANIONS  2  1 68E  -  03 

PERCENT  DIFFERENCE  ■  -  43  (CATIONS  -  AN1 ONS > / < AN  I ONS  ♦ 

NONCARBONATE  ALKALINITY  «  0 

IONIC  STRENGTH  *  3  026E - 03 


CATIONS) 


Table  10 

I'CNTBQ  Output 
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is  selected  is  dependent  upon  the  MISP  file  configura¬ 
tion.  The  parameter  which  controls  the  EXAMS  output  is 
the  number  of  times  EXAMS  uses  MINTEQ  to  update  steady 
state  metal  concentrations.  EXAMS  utilizes  MINTEQ  three 
times  to  update  steady-state  metal  concentrations  in  the 
Naugatuck  River  test  case.  The  selection  of  this 
parameter  is  based  on  trial  and  error;  the  programmer 
examines  the  output  obtained  while  varying  the  input  and 
determines  when  system  steady  state  is  achieved.  System 
steady  state  is  defined  as  an  equilibrium  between  the 
partioning  achieved  with  EXAMS  and  MINTEQ. 

The  EXAMS  output  file  received  after  each  steady 
state  calculation  will  be  divided  into  three  sections. 

The  first  section  consists  of  13  tables,  which  were 
obtained  from  the  initial  EXAMS  calculation.  The  initial 
EXAMS  calculation  assumes  the  pollutant  is  in  the  aqueous 
form.  Tables  2-8  (all  shown  in  Table  11)  are  echoes  of 
the  EXAMS  input  data.  Table  1  is  not  printed  because  the 
pollutant  in  question  is  a  heavy  metal  and  not  an  or¬ 
ganic.  If  the  pollutant  had  been  an  organic,  Table  1 
would  have  been  an  echo  of  the  organic  input  data. 

The  next  section  includes  Tables  9-14  (all 
shown  in  Table  12).  These  five  tables  are  the  transport 
profile,  the  kinetic  profile,  the  canonical  profile, 
toxicant  loadings  by  system  element,  the  distribution  of 
the  chemical  at  steady  state  (the  water  column  and 
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sediments),  the  average,  maximum  and  minimum  concentra¬ 
tions  at  steady  state.  Table  10,  the  kinetic  profile  is 
for  organic  toxicants  only.  These  five  tables  are 
repeated  the  number  of  times  EXAMS  uses  MINTEQ  to  update 
steady  state  metal  concentrations.  For  the  Naugatuck 
River  test  case,  this  occurred  three  times.  The  most 
important  of  the  three  outputs  is  the  last  repetition, 
because  it  is  the  output  obtained  when  system  steady 
state  has  been  achieved  for  the  water  column  and  sedi¬ 
ments. 

The  third  section  of  output  is  encompassed  in 
Table  15  -  17,  the  analysis  of  steady  state  fate,  the 
simulation  of  system  response  after  load  ceases,  and  the 
exposure  analysis  summary.  These  three  tables  contain  the 
persistence  calculation  performed  by  EXAMS.  These  tables 
are  presented  in  Table  13.  The  final  portion  of  the 
program  output  is  the  MISP  output. 

3.5  Output  {MISPj 

This  section  as  previously  stated  provides  a 
summary  of  the  EXAMS-MINTEQ  interactions.  These  interac¬ 
tions  include  steady  state  calculations  made  to  achieve 
system  steady  state,  and  calculations  made  during  the 
persistence  to  update  metal  fractions.  For  the  Naugatuck 
River  test  case,  steady  state  concentrations  were  updated 
three  times  to  obtain  system  equilibrium,  and  updated 
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Table  14 


MISP  Output 


UPDATED  CONCENTRATIONS  FROM  MINTED  (MG/L) 
UPDATED  CONCENTRATIONS  FROM  MINTEQ  (MG/L) 


COMPARTMENT 

AQUEOUS 

ADSORBED 

PRECIPITATED 

COMPARTMENT 

AQUEOUS 

ADSORBED 

PRECIPITATED 

1 

9. 024E-03 

2. 478E-03 

0.0 

2 

9. 532E-03 

1. 379E+02 

0.0 

3 

1.341E-02 

1. 519E-03 

0.0 

4 

1 . 836E-02 

5. 897E+0 1 

0.0 

5 

1. 684E-02 

7. 880E-03 

0.0 

6 

1 . 697E-02 

7.435E+01 

0.  0 

7 

1.404E-02 

2 . 065E-03 

0.0 

8 

5.420E-03 

1. 798E+02 

0.0 

9 

1. 061E-02 

1. 970E-03 

0.0 

10 

6. 094E-03 

1. 500E+02 

0.0 

CONCENTRATION 

IS  COMPUTED 

BY  EXAMS  IN 

MG/L 

CONCENTRATIONS  COMPUTED 

BY  EXAMS  IN 

MG/L 

'  V  1L*.  k. 
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CHAPTER  4 


MODEL  RESULTS 

4.1  Calibration  Parameters 

Metal  transport  in  the  Naugatuck  River  is  cur¬ 
rently  being  investigated  to  provide  important  data  for 
waste  allocation  decisions.  Included  within  this  inves¬ 
tigation  is  the  calibration  of  MEXAMS  to  a  compartmental¬ 
ized  river  reach.  The  reach  being  calibrated  (Figure  3j 
was  originally  composed  of  three  water  column  compart¬ 
ments  and  three  benthic  compartments.  However,  the  final 
model  configuration  consisted  of  a  ten  compartment  model; 
five  benthic  compartments  and  five  water  column  compart¬ 
ments.  A  diagram  of  the  final  configuration  will  be 
presented  herein.  This  calibration  consisted  of  adjusting 
certain  model  inputs  to  yield  results  that  were  the  same 
as  the  field  data  obtained  on  August  16-18,  1983. 

The  first  step  in  the  calibration  was  to  deter¬ 
mine  which  model  inputs  have  the  least  reliability.  These 
are  the  parameters  that  will  be  adjusted  to  perform  a 
sensitivity  analysis.  Any  sensitivity  analysis  is  system 
dependent  because  certain  input  parameters  may  be  more  or 
less  sensitive  depending  on  the  system  being  modeled. 


Therefore,  the  first  decision  to  be  made  was  to  select  a 
system  conf  iguration.  There  are  many  possible  configura¬ 
tions  due  to  the  possible  advective  and  dispersive  path¬ 
ways.  The  two  initial  configurations  chosen  for  this  test 
case  are  presented  in  Figure  5.  The  main  difference 
between  configuration  one  and  configuration  two  is  the 
advective  pathways.  Configuration  one  allows  flow  to  exit 
the  system  through  the  sediments  via  compartment  six  and 
has  no  advective  flow  through  compartment  four,  con¬ 
figuration  two  is  just  the  opposite.  Originally,  both  of 
these  configurations  had  the  entire  groundflow  routed 
through  the  benthic  compartments  (2  &  6).  This  caused  a 
considerable  problem  because  the  sediment  residence  time 
within  the  benthic  compartments  was  so  low  that  the 
priority  pollutant  was  being  flushed  out  of  the  system. 
This  problem  was  corrected  by  routing  a  portion  of  the 
groundwater  into  the  water  column  compartments  (1  &  5). 
This  increased  the  sediment  residence  time  and  thus, 
increased  the  concentration  of  adsorbed  pollutant  within 
the  bed  sediments.  Figure  6  displays  the  effect  of  rout¬ 
ing  various  percentages  of  groundwater  through  the  ben¬ 
thic  compartments  for  both  configurations.  Compartments 
two  and  four  were  not  shown  in  configuration  one  because 
the  residence  time  changed  very  little.  This  slight 
change  is  a  result  of  no  flow  being  advected  through 
these  compartments.  On  the  other  hand,  compartment  six  is 
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CONFIGURATION  I 


J  LU 


CONFIGURATION  t 
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Figure  5. 
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not  shown  for  configuration  two  because  it  also  had  very 
little  change.  Again,  this  is  a  result  of  no  flow  being 
advected  through  compartment  six.  At  this  point  a  deci¬ 
sion  was  made  to  perform  the  remaining  sensitivity 
analysis  and  calibration  on  configuration  two. 

The  effect  of  the  percent  groundwater  flow  into 
the  benthic  compartments  for  configuration  two  allowed 
adjustment  of  the  particulate  concentration  in  the  ben¬ 
thic  regions  (Figure  7).  The  next  parameter  adjusted  for 
configuration  two  was  the  dispersion  coefficient.  This 
parameter  was  adjusted  from  5.0Q*E-05  to  5.00*E-08 
(meters  square/hr)  to  obtain  the  relationship  between 
dispersion  coefficient  and  adsorbed  concentration  in  the 
benthic  compartments  (Figure  8).  The  above  range  for 
dispersion  coefficients  did  not  result  from  a  calculation 
because  the  slope  of  the  river  bed,  the  magnitude  of  the 
diffusive  exchange  coefficients  for  the  water  column  and 
the  interstitial  water,  and  the  characteristic  length  of 
the  gradient  were  not  known.  However,  a  range  of  possible 
dispersion  coefficients  was  obtained  from  Analysis  of 
Fate  of  Chemicals  in  Receiving  Waters  -  Phase  I  (1981, 
DiToro  and  St.  John;. 

The  trend  shown  by  varying  the  dispersion  coeffi¬ 
cient  resulted  in  the  sensitivity  analysis  of  a  final 
parameter.  This  parameter  was  the  and  originally  it 
was  not  to  be  a  calibration  variable.  The  reason  behind 
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DISPERSION  COEFFICIENT  VS  CONCENTRATION 
IN  BENTHIC  COMPARTMENTS 


Figure  8. 


the  change  stemmed  from  the  inability  to  achieve  the 
required  adsorbed  concentration  within  the  benthic  com¬ 
partments.  The  decision  to  increase  the  «d  values  from 
their  field-determined  values  was  reasonable  considering 
the  accuracy  of  benthic  Kd  determinations  made  in  the 
field.  It  is  very  difficult  to  obtain  a  sediment  sample 
without  entrainment  of  water  from  the  water  column.  The 
aqueous  concentration  of  contaminant  within  the  water 
column  is  much  lower  than  that  within  the  bed  sediments, 
and  the  result  would  be  a  dilution.  This  dilution  wou7d 
result  in  a  value  that  was  too  low,  and  for  this 
reason  the  parameter  has  been  changed  to  achieve 
higher  particulate  concentrations  within  the  bed  sedi¬ 
ments.  The  adjustment  of  the  value  and  its  affect  on 
the  benthic  compartment  adsorbed  concentrations  are  shown 
in  Figure  9. 

Another  problem  was  encountered  after  the 
calibration  of  configuration  two  which  resulted  in  a 
final  adjustment  of  the  advective  pathways.  The  problem 
was  that  the  total  concentration  in  compartment  five  was 
too  high,  and  to  reduce  this  an  advective  pathway  from 
the  benthic  compartment  (compartment  six)  was  es¬ 
tablished.  The  advective  flow  proportions  from  compart¬ 
ment  six  to  compartment  five  were  varied  in  effort  to 
reduce  the  total  concentration  in  compartment  five.  The 
results  obtained  from  varying  the  advective  flow  are 
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presented  in  Figure  10. 

The  required  concentrations  for  compartment  five 
could  not  be  obtained  by  varying  the  advective  flow.  As  a 
result,  a  ten  compartment  configuration  was  developed. 
Since  the  model  establishes  instantaneous  equilibrium, 
the  6.7  mile  compartment  was  to  large  and  the  concentra¬ 
tion  could  not  be  lowered  by  adjusting  any  of  the 
calibration  parameters.  The  final  ten  compartment  con¬ 
figuration  is  presented  in  Figure  11.  The  final  con¬ 
figuration  was  calibrated  by  adjusting  three  variables; 
the  percent  ground  water  flow  into  the  benthic  compart¬ 
ments  v8  &  10;,  the  advective  flow  between  compartments 
(8  to  7,  and  10  to  9),  and  the  Kd  value  for  the  benthic 
compartments. 

4.2  Model  Results 

After  a  good  understanding  was  achieved  for  the 
sensitivity  of  the  calibration  parameters,  the  calibra¬ 
tion  continued  for  the  ten  compartment  model.  The  two 
models  being  compared  are  slightly  different,  because  the 
incremental  flow  model  utilizes  point  source  loads  as 
incremental  loadings  to  the  an  existing  concentration. 
This  new  concentration  is  then  transported  downstream.  On 
the  other  hand,  the  MEXAMS  model  establishes  instan¬ 
taneous  equilibrium  with  the  point  source  load,  and  as  a 
result,  the  model  outputs  are  mean  compartment  outputs 
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%  ADVECTED  FLOW 


and  not  outputs  at  various  points  within  the  r 
Thus,  the  MEXAMS  model  will  not  predict  the  cc 
where  a  point  load  enters  iBogue  Rd.)>  and  the 
tion  measured  will  be  higher  than  predicted.  1 
calibration  was  compared  to  an  incremental  flc 
Curtis,  State  of  Connecticut,  Dept,  of  Environ 
Protection/  and  laboratory  results  obtained  fr 
collected  at  various  sites  along  the  river.  Th 
compartment  model  calibration  to  particulate  c 
tions  within  the  benthic  compartments  was  very 
Tne  model  calibration  to  aqueous  concentratior 
the  water  column  was  not  as  accurate  because  c 
five  concentrat ion  were  too  high.  As  a  result, 
compartment  was  further  divided  into  smaller  c 
ments.  The  six  compartment  model  results  are  p 
Figure  12  and  Table  15. 

This  model  was  improved  by  dividing  tl 
river  segment  into  smaller  compartments.  The  c 
talization  of  compartments  five  and  six  into  ; 
compartments  would  allow  for  a  greater  control 
model  calibration.  The  third  segment  has  been 
figured  by  splitting  compartments  five  and  si) 
The  calibration  results  obtained  from  the  fine 
figuration  are  presented  in  Figure  13  and  Tab' 
respectively.  As  stated  previously,  three  par; 
adjusted  to  achieve  a  calibrated  model.  Table 
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Table  17 

Comparative  Kd  Data 

Compartment  i  Original  log  Kd  Calibrated  log  Kd 


1 

0.30 

0.30 

2 

3.72 

5.02 

3 

-0.25 

-0.25 

4 

2.95 

4.20 

5 

-0.05 

-0.05 

b 

2.67 

4.42 

7 

-0.05 

-0.05 

8 

2.67 

5.30 

9 

-0.05 

-0.05 

10 

2.67 

5.17 

the  original  log  s  and  the  final  log  K^'s  used  to 
obtain  the  ten  compartment  calibration.  The  advective 
flow  pa*hs  for  the  new  compartmental iz ations  were 
presented  in  Figure  11,  and  the  percent  ground  water 
routed  into  benthic  compartments  eight  and  ten  were 
changed  from  .01%  to  .1%  and  .13%,  respectively.  The 
final  calibration  for  the  ten  compartment  model  was  very 
accurate,  and  the  concentrations  predicted  by  the  model 
were  well  within  the  error  bars  of  the  field  data. 


4.3  Future  Loadings 


Since  the  purpose  of  this  report  is  to  illustrate 
the  usefulness  of  the  modeling  system,  the  effect  of  an 
increase  in  metal  loading  to  compartment  one  was  simu¬ 
lated.  This  increase  could  be  an  additional  discharger  at 
some  point  in  the  future,  or  an  increase  in  production  by 
an  existing  discharger.  The  existing  contaminant  loading 
is  1.284  E - 02  Kd/hr  and  this  was  increased  to  1.284  E-01 
and  1.284  E+00.  The  effect  of  this  increased  metal  load¬ 
ing  on  both  water  column  and  benthic  compartments  is 
displayed  in  Figure  14.  This  model,  however,  has  some 
limitations,  and  they  will  be  discussed  next,  along  with 
recommendations  for  solving  these  drawbacks. 


4.4  Limitations  and  Recommendations 


The  first  limitation  encountered  dealt  with  the 


COMPARTMENT  COMPARTMENT 
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ADSORPTION  IN  COMPARTMENT  I. 
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Other  limitations  include  the  ability  to  only  run 
the  program  with  one  metal.  This  allows  no  competition 
and  is  unrealistic  for  an  aquatic  ecosystem..  Also,  on 
the  subject  of  competition,  when  using  the  «d  "activity" 
approach,  unlimited  surface  sites  are  assumed  which  is 
unrealistic.  It  is  recommended  that  an  input  be  es¬ 
tablished  that  would  limit  surface  sites  and  make  the 
model  more  realistic. 


CHAPTER  5 


CONCLUSIONS 

The  behavior  of  total,  dissolved  and  particulate 
(adsorbed  copper  within  an  8.3  mile  reach  of  the 
Naugatuck  River  (East  Albert  St.  to  Thomaston  Bridge 
Abutment)  was  modeled  using  the  Metals  Exposure  Assess¬ 
ment  Modeling  System  ^MEXAMS).  Data  provided  by  the 
Department  of  en v i ronmenta 1  Protection,  State  of  Connec¬ 
ticut  concerning  mass  metal  loading,  system  physical  and 
chemical  characteristics  was  used  to  calibrate  the  MEXAMS 
computer  model  for  this  test  case. 

The  most  important  aspect  of  this  model  calibra¬ 
tion  was  the  compartmental ization  (configuration)  of  the 
river  reach.  The  model  was  initially  calibrated  for  a  six 
compartment  configuration.  However,  the  six  compartment 
configuration  did  not  predict  adequately  (within  10% 
error)  the  total  and  dissolved  concentrations  in  the 
furthermost  downstream  water  column  compartment  (compart¬ 
ment  five).  The  inability  for  the  six  compartment  model 
to  adequately  predict  the  target  (field)  concentrations 
in  compartment  five  was  due  to  the  length  of  the  compart¬ 
ment  (6.7  miles),  and  the  fact  that  the  model  assumes 
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complete  mixing  and  creates  instantaneous  equilibrium 
within  each  compartment.  As  a  result,  a  ten  compartment 
configuration  was  utilized  to  obtain  the  final  calibra¬ 
tion.  This  configuration  was  obtained  by  dividing  com¬ 
partments  five  (water  column)  and  six  (adjacent  benthic) 
into  three  equal  compartments.  This  division  was  per¬ 
formed  to  achieve  a  greater  compartment  sensitivity  to 
obtain  the  target  (field  measured)  concentrations. 

This  final  conf iguration  abased  upon  the  data 
provided)  accurately  predicted  total,  dissolved  and 
particulate  concentrations  for  all  water  column  and 
benthic  compartments.  The  major  parameters  adjusted  to 
obtain  this  final  calibration  were  the  percent 
groundwater  routed  through  the  benthic  compartments,  the 
percent  advected  flow  from  benthic  to  water  column  com¬ 
partments,  and  the  "activity"  Kd’s  were  kept  at  field 
determined  values.  However,  the  calibration  of  this  model 
is  only  the  initial  phase  of  accurately  modeling  this 
river  reach,  and  model  validation  to  another  data  set  is 
required. 

The  validation  process  could  also  include  a 
reduction  of  the  following  model  limitations  that  were 
discovered  during  the  calibration  process;  the  need  for 
additional  organic  complexes  within  the  thermodynamic 
data  base,  the  ability  for  the  model  to  simulate  more 
than  one  priority  pollutant  metal  simultaneously,  and  the 
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ability  for  the  user  to  adjust  the  concentration  to  be 
obtained  during  the  system  self-purification  (persistence 
calculation).  The  validation  of  this  model  could  result 
in  possible  future  applications  and  provide  additional 
assistance  to  the  Department  of  Environmental  Protection, 
State  of  Connecticut  in  determining  waste  discharge 
allocations  for  present  and  future  dischargers  within 
this  reach  of  the  Naugatuck  River. 

In  conclusion,  the  results  of  the  study  suggest 
that  MEXAMS  is  a  potentially  powerful  management  tool 
that  could  be  utilized  for  a  variety  of  environmental 
problems.  Future  activity  using  the  model  should  continue 
the  application  to  actual  field  situations  so  that  con¬ 
tinued  refinement  of  the  model  and  appropriate  modifica¬ 
tions  are  made  to  promote  its  use. 
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APPENDIX  A 


DESCRIPTION  OF  THE  MINTEQ  INPUT  FILE 
(Felmy  et  al . ,  1 983 j 

This  appendix  describes  how  to  setup  the  sample 
description  file  for  running  the  VAX  and  PDP  11/70  ver¬ 
sions  of  MINTEQ.  A  detailed  description  of  the  options 
available  to  the  user  accompanies  the  file  descriptions. 

The  file  is  broken  into  three  sections.  Each 
section  is  separated  by  a  blank  line.  There  should  be 
three  and  only  three  blank  lines  in  the  sample  file.  The 
three  sections  will  be  called  Basic  Input  Data,  Type 
Modifications  and  Insertion  of  Species.  The  format  field 
designations  for  each  line  are  given  in  parenthesis 
following  the  line  designation. 

Section  1  -  Basic  Input  Data 

Line  1 iA80j.  Alphanumeric  description  of  the  sample. 

Line  2  (AbO j . Alphanumeric  description  of  the  sample. 

Line  3 . _ Water  temperature,  analytical  units  and  fixed 

ionic  strength  designation 

The  water  temperature  in  degrees  Celcius  is  entered 
followed  by  the  analytical  units.  In  the  VAX  free-format 
version,  the  temperature  and  units  must  be  separated  by  a 
space  and  the  units  enclosed  in  quote  marks.  The  PDP 
11//0  version  follows  the  designated  format  fields.  The 
analytical  unit  designations  can  be  either  'PPM',  ’ M G / 1 '  , 
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'MOL',  or  'MEQ/L'.  These  designations  stand  for  parts  per 
million,  milligrams  per  liter,  molality  and  mill i - 
equivalents  per  liter,  respectively.  The  fixed  ionic 
strength  is  entered  following  the  analytical  units.  In 
the  VAX  version  there  must  be  a  space  between  the 
analytical  units  and  the  ionic  strength.  If  the  ionic 
strength  is  not  fixed,  enter  0.0  following  the  analytical 
units. 

Line  4L8  i  1 1 , 1 X ) j .  This  line  contains  the  run  specific 
user  options. 

Option  1,  The  inorganic  carbon  input  option. 

(0)  =  Total  inorganic  carbon. 

(l)  =  Total  alkalinity. 

Option  2.  The  Debug  print  option.  This  is  generally 
used  only  when  modifications  have  been  made  to  the 
code  and  the  values  of  certain  arrays  must  be 
checked. 

(0)  =  No  Debug  printout. 

(1)  =  Prints  the  value  of  the  IDX,  C,  IDY, 

GAMMA,  C  and  Y  arrays. 

(3)  =  Prints  the  NNN,  NN(1),  NN(2)  and  the 

IDY,  GK,  GK 1  arrays  in  subroutine  KCORR. 

(4)  =  Prints  the  IDX,  C,  IDY,  L  and  GK 

arrays  every  iteration  in  subroutine  SOLID 
and  the  saturation  index  for  the  solid 
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added  to  the  phase  assemblage  in  sub¬ 
routine  SOLIOX. 

Option  3.  The  charge  balance  option. 

(0)  =  Terminates  execution  of  the  program  if  the 
initial  charge  in  balance  is  greater  than 
30  percent. 

(1)  =  Does  not  terminate  execution  regardless 
of  the  initial  computed  charge  imbalance. 

Option  4.  Considered  solids  and  print  option. 

(0)  =  Do  not  allow  all  of  the  solids  in  the 

data  base  to  precipitate  or  dissolve.  The 
only  solids  considered  will  be  those 
entered  in  the  next  input  section  under 
type  changes.  Print  the  problem  results 
after  the  initial  aqueous  speciation  plus 
solids  problem  is  solved  and  after  all 
Type  V  solids  are  either  in  equilibrium  or 
under  saturated. 

(lj  =  Allow  all  solids  in  the  data  base  to 

precipitate  if  they  become  over  saturated. 

That  is,  designate  all  solids  in  the  data 
base  as  Type  V.  Print  the  problem  results 
only  after  the  entire  problem  has  been 
solved. 


(2)  =  Consider  all  solids  in  the  data  base. 


Print  the  problem  results  after  the  ini¬ 
tial  user  specified  problem  has  been 
solved  and  again  after  all  Type  V  solids 
are  undersaturated  or  in  equilibrium. 

(3)  =  Consider  all  solids  in  the  data  base. 

Print  the  problem  results  following  the 
selection  of  every  solid  and  after  all  of 
the  solids  are  in  equilibrium  or  under¬ 
saturated. 

Option  5,  The  total  number  of  iterations  option. 

(0)  =  Allow  40  iterations. 

(1)  =  Allow  10  iterations. 

(2)  =  Allow  100  iterations. 

(3)  =  Allow  200  iterations. 

The  100  or  200  iteration  options  should  be  selected 
if  a  large  number  of  solids  have  been  designated 
as  considered  solids.  The  ten  iteration  option  is 
only  useful  when  debugging  the  program. 

Option  6.  The  pH  variation  option. 

(0)  =  Do  not  allow  the  pH  to  automatically 

vary  during  precipitation/dissolution  of 
sol  id  phases. 

(1)  =  Allow  the  pH  to  vary  during  precipita¬ 


tion/dissolution  of  solid  phases.  Not  for 
this  option  to  work  all  solid  phases  must 
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be  declared  Type  V  in  the  initial  input. 

Option  7.  The  fixed  ionic  strength  option. 

(0)  =  Allow  MINTEQ  to  compute  the  ionic 
strength. 

(1)  =  Fix  the  ionic  strength  at  the  value 
designated  on  line  3. 

Option  8.  The  numerical  method  option. 

(0)  =  Use  only  Newton-Raphson  iteration. 

(1)  =  Use  a  combination  of  Newton-Raphson 
and  a  modified  line  search.  This 
option  should  only  be  used  after 
consulting  the  technical  manual. 

Option  9.  Output  Option  (available  only  on  VAXj 
(0)  =  Do  not  print  on  initial  listing  of  the 
thermodynamic  data. 

(1)  =  Print  a  listing  of  the  thermodynamic  data. 

Line  5  Li 1 , IX ,4 (F6 . 2 ,  IX )  |.  This  line  is  for 
input  of  adsorption  parameters  other  than  the  mass  total 
and  activity  guesses.  The  first  parameter  on  this  line  is 
the  adsorption  model  being  used  (II  field).  The  options 
are: 

0  -  No  Adsorption 

1  -  Activity  Kd,  langmiur  isotherm 
Freundlich  isotherm 


2  -  Constant  Capacitance  Double  Layer  Model 

3  -  Triple  Layer  Site  Binding  Model. 

The  next  four  inputs  on  this  line  are: 

-  Solid  Concentrat ion  (g/1) 

-  Specific  Surface  Area  im2/g) 

-  Inner  Layer  Capacitance  (F/m2) 

-  Outer  Layer  Capacitance  (F/m2) 

input  in  this  order.  None  of  the  last  four  inputs  are 
required  for  the  activity  Kd.  The  constant  capacitance 
model  does  not  require  an  outer  layer  capacitance.  The 
triple  layer  model  requires  all  four  data  inputs. 

Line  6+  ( 1 7  ,  IX ,  E9 . 3 , 1 X  ,  F 6 . 2 ) .  Component  input  lines. 
There  are  as  many  of  these  lines  as  there  are  components. 
A  blank  line  must  follow  the  last  component. 

The  first  specification  (17)  is  the  component  ID 
number.  Component  ID  numbers  are  given  in  Table  3.  The 
second  specification  (E9.3  field)  is  the  total  analytical 
mass  in  the  units  designated  on  line  3.  The  only  excep¬ 
tion  occurs  when  the  constant  capacitance  of  triple  layer 
models  are  used.  In  such  cases  the  analytical  mass  for 
the  surface  sites  (S0H1  or  S0H2  components)  must  be 
specified  in  sites  per  gram.  The  third  specification  on 
this  line  is  a  guess  at  the  log  of  the  component  ac¬ 
tivity.  If  you  leave  this  field  blank  the  initial  es¬ 
timate  of  the  activity  will  be  the  analytical  molality 
divided  by  100.  Remember  you  need  one  line  for  every 


component  and  a  blank  line  after  the  last  component. 


Section  2  -  Type  Modification 

This  section  is  for  changing  the  default  species 
designations.  The  default  type  specification  have  been 
described  in  Section  5. 

The  first  line  of  this  section  contains  the  first 
species  type  designation  and  the  number  of  species  of 
this  type  Format  (13,  IX,  13;.  The  type  designations  can 
range  from  two  to  six.  Then  for  each  species  of  the 
entered  type,  a  line  is  included  (17,  IX,  E9.3,  IX,  F6.2; 
designating  information  for  that  species.  The  first  field 
specification  is  for  the  species  ID  number.  Species  ID 
numbers  for  components  (Type  I)  are  given  in  Table  3.  All 
other  species  ID  numbers  can  be  found  in  the  listing  of 
Thermodynamic  Data  given  in  the  MINTEQ  Technical  Report. 
The  next  two  specifications  are  for  the  new  log  K  and 
enthalpy  of  reaction.  These  are  both  optional  and  if  not 
included  the  default  values  in  the  thermodynamic  data 
base  are  used.  A  blank  line  also  ends  this  section.  In 
the  case  of  Type  IV  species  with  an  initial  mass  total 
(in  moles/1)  there  is  an  additional  input  field  for  the 
initial  mass.  This  input  field  follows  the  input  of  the 
enthalpy  of  reaction  resulting  in  four  inputs  on  one 
line.  The  input  format  for  Type  IV  species  is  (17,  IX, 
E9.3,  IX.  F6.2,  IX,  E10.3).  For  examples  of  necessary 
species  modifications  to  solve  specific  problems  see 


Section  5. 


Section  3  -  Insertion  of  Species  Not  in  the  Data  Base 
The  firrt  line  is  for  designating  the  species  type 
and  the  number  of  new  species  of  this  type  (Format  13, 

IX,  13).  The  species  type  can  only  range  from  two  to  six. 
The  next  lines  contain  the  data  for  the  new  species  of 
the  specified  type.  There  are  three  lines  for  each 
species.  The  first  line  of  species  data  is  in  Format  (17, 
IX,  A12 ,  2F10.4,  2F8.3,  3F5.2,  F9.4).  The  format  fields 
correspond  to  the  following  data:  ID  number,  name,  en¬ 
thalpy  of  reaction,  log  K,  minimum  log  K,  maximum  log  K, 
charge,  Debye-Huckel  A  parameter,  Debye-Huckel  B 
parameter  and  molecular  weight  respectively.  Only  the  ID 
number  and  log  K  are  absolutely  essential.  For  a  descrip¬ 
tion  of  the  Debye-Huckel  parameters  see  the  MINTEQ  Tech¬ 
nical  Report.  The  second  line  of  data  is  in  Format  LFB.2, 
IX,  II,  IX,  6 (F7 . 3 ,  IX,  13,  IX)].  The  format  fields 
correspond  to:  carbonate  alkalinity  factor,  number  of 
components  in  the  reaction  and  the  stoichiometry  and  ID 
number  for  up  to  six  components.  The  carbonate  alkalinity 
factor  is  only  useful  if  the  input  inorganic  carbon  is  an 
alkalinity  value.  A  description  of  the  carbonate 
alkalinity  factor  is  given  in  the  MINTEQ  Technical 
Report.  The  third  line  is  in  Format  L3(F7.3,  IX,  13,  lX;j 
and  is  merely  a  continuation  of  the  component  entries  for 
the  inserted  species.  In  the  VAX  version,  the  third  line 


98 


is  in  format  [3(F7.3,  IX,  13,  IX),  II,  3(F7.3,  IX,  13, 
1X)j.  The  format  field  beginning  with  II  is  for  inserting 
species  information  into  the  "B"  matrix.  Such  information 
is  only  useful  when  the  component  stoichiometries  in  the 
mass  action  expressions  are  different  from  the 
stoichiometries  in  the  mass  balance  equation.  The  latter 
format  fields  correspond  to  the  number  of  components  in 
the  "B"  matrix,  stoichiometry  of  the  component  and  com¬ 
ponent  ID  number.  This  section  is  also  terminated  with  a 
blank  line. 


APPENDIX  B 


DESCRIPTION  OF  THE  EXAMS  INPUT  FILE 
{Burns  et  al . ,  1 981 , 


The  narrative  associated  with  each  database  references 
the  format  codes  listed  below. 


Code 

1 

2 

3 


4 

5 


Format 

60  characters,  left  justified. 
Format:  { 6 0 A 1 ) 

Columns  1-5  right  justified. 

Format:  (2015) 

Column  1,  one  character/record.  The 
maximum  number  of  records  is  equal 
to  the  number  of  compartments. 
Format:  (Al) 

Columns  1-10.  Format:  (FlO.Oj 
Each  record  contains  8  10-column 
fields,  maximum  2  records  per 
variable.  Format:  (8F10.0) 


Record 


Columns  Compartment 

1-10  1 

11-20  2 


Record  2 
(if  needed) 


1-10 

11-20 


9 

10 


6  Five  record  layout  containing  39 

values.  Four  records  contain  eight 
10-column  fields:  the  fifth  contains 
seven.  Format:  (8F20.0) 


Columns  Value 

Record  1  1-10  1 

11-20  2 


71-80  8 

Record  2  1-10  9 


Record  5 


1-10 


33 


61-70  39 


The  number  of  records  is  equal  to 
NSPEC1  or  NSPEC2  divided  by  16.  Each 
record  contains  16  five-column 
fields.  For  JFRAOG  and  ITOADG  use 
NSPEC1  and  for  JTURBG  and  ITURBG  use 
NSPEC2.  Each  field  is  right  justi¬ 
fied.  Format:  (2015) 


Columns  <Par> 

Record  1  1-5  1 

6-10  2 


76-80  16 

Record  2  1-5  17 


•  • 

Record  n  41-50  29 

51-60  30' 

Format:  (5(11 ,1X) ) 

Column  Flag 

1  1 

3  2 

5  3 

7  4 

9  5 

Format:  (8F10.0) 

Columns  Value 


Format:  (8F10.0) 


Columns 

1-10 

11-20 

21-30 

31-40 


Value 

1 

2 

3 

4 


Format:  (8F10.0) 

Columns 

1-10 

11-20 


Value 

1 

2 


The  number  of  records  is  equal  to 
NSPEC1  or  NSPEC2  divided  by  eight. 
Each  record  contains  eight  10-column 
fields.  For  ADVPRG  use  NSPEC1  and 
for  XSTURG,  CHARLG,  and  DSPG  use 
NSPEC2.  Format: (8F10.0) 


Columns 

Record  1  1*10 

11-20 


Value 


1 

2 


Recori  2 


71-80 


41-50  29 

51-60  30 


Record  n 


Creating  or  Changing  the  User  Run  Information 


This  sequential  file  contains  chemical  name,  ecosystem, 
name,  and  loads. 


Record 

1: 

(CHEM1 ) 

Chemical  name  -  three  characters 
maximum.  Format: (A3) 

Record 

2: 

(EC01 ) 

Ecosystem  name  -  three  characters 
maximum.  Format:  (A3) 

Record 

3: 

(NLOAD) 

Number  of  Loads.  Format:  (15) 

Record 

4,. 

. . ,NL0AD+3:(I,  STRLDG,  NPSLDG,  PCPLDG,  DRFLDG , 

IFLLOGJ  Each  of  NLOAD  records 
contain  the  load  number  and  the 
I th  STRLDG,  NPSLDG ,  PCPLDG, 

ORFLOG,  IFLLDG.  Format 
(I5.5F10.0; 

STRLDG  -  Stream  loadings. 

NPSLDG  -  Non-point  source  loadings. 
PCPLDG  -  Precipitation  loadings. 
ORFLDG  -  Drift  loadings 
IFLLDG  -  Interflow  loadings. 


Record  4 


Columns 

1-5 

6-15 

16-25 

26-35 

36-45 

46-55 


Field 

I 

STRLDG(I) 
NPSLDG ( I ; 
PCPLDG( I) 
DRFLDG(I) 
IFLLDG(I) 
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Creating  or  Changing  the  Environment  Database 


The  compartment  parameter  data  are  located  in  the  follow¬ 
ing  sequential  record  layout.  (The  format  code  list 
provides  detailed  format  specifications.) 


Record : 

( SYSTYPj 

Ecosystem  name,  maximum  60  characters 
Format:  1 

Record: 

(KOUNT) 

Number  of  compartments,  maximum  number 
of  compartments  is  NPX,  which  is 
fixed  for  a  particular  installation. 
Format:  2 

Record: 

(TYPEEj 

The  next  KOUNT  records  each  contain  a 
one  character  compartment  code, 

KOUNT  the  number  of  compartments. 

Format:  3 

Record: 

(LATG) 

Ecosystem  latitude,  columns  1-10, 
degrees  and  seconds  as  whole  number 
and  decimal  fraction,  e.g.  10.15. 
Format:  4 

Record: 

(WINDG; 

Wind  velocity  per  compartment,  eight 
per  record.  Format:  5 

Record: 

( VOLG) 

Total  environmental  volume  per  compart¬ 
ment.  Format:  5 

Record : 

(AREAG) 

Area  of  each  ecosystem  compartment. 
Format:  5 

Record: 

l  DEPTHG ) 

Average  depth  of  each  compartment. 
Format:  5 

Record : 

( STDLOG) 

Stream  flow  entering  ecosystem  compart¬ 
ments.  Format:  5 

Record: 

(STSEDG) 

Stream-born  sediment  load  entering 
ecosystem  compartments.  Format:  5 

Record: 

( NPSFLGj 

Non-point-source  water  flow  entering 
ecosystem  compartments.  Format:  5 

Record: 

(NPSEDG) 

Non-point-source  sediment  loads  enter¬ 
ing  ecosystem  compartments.  Format:  5 

Record : 

(INTFLG) 

Interflow  (groundwater  seepage)  enter- 

ing  each  ecosystem  compartment. 


Format:  5 


Record:  ^RAlNGj  Average  rainfall  in  geographic  area  of 
system.  Columns  1-10.  Format:  4 

Record:  (C10UDG)  Average  cloudiness  in  tenths  of  full 
sky  cover.  Columns  1-10.  Format:  4 

Record:  (WLAMG)  Temporally  averaged  spectral  irradiance 
immediately  below  the  water  surface. 
This  is  five  record  layout  containing 
39  values,  ten  columns  per  value, 
eight  values  per  record.  Format:  6 

Record:  (DFACG)  Distribution  function  (optical  path) 
for  each  compartment.  Format:  5 

Record:  (EVAPG)  Evaporative  water  losses  from  ecosystem 
compartments.  Format:  5 

Record:  (SDCHRG)  For  water  column  compartments  the 

suspended  sediment  concentration;  for 
benthic  sediment  compartments  the  bulk 
density  of  the  bottom  sediment. 

Format:  5 

Record:  (PCTWAG)  Percent  water  in  bottom  sediments  of 
benthic  compartments.  Elements  corres¬ 
ponding  to  water  column  compartments 
are  not  used  (dummy  values).  Format:  5 

Record:  *.NSPEC1J  Number  of  active  advective  transport 
pathways.  Format:  2 

Record:  (JFRADGj  Source  compartment  [ J )  for  advective 
flow.  Format:  7 

Record:  (ITOADG)  Receiving  compartment  (ij  for  advective 
flow.  Format:  7 

Record:  (ADVPRGj  Proportion  of  total  advective  flow  from 
compartment  J  that  flows  to  compart¬ 
ment  I.  Format:  12 

Record:  (NSPEC2)  Number  of  active  dispersive  transport 
pathways.  Format:  2 

Record:  (JTURBG)  Source  compartment  for  dispersive  flow. 
Format:  7 

Record:  (ITURBGj  Receiving  compartment  for  dispersive 
flow.  Format:  7 


Record:  (XSTURG)  Cross-sectional  areas  for  dispersive 
exchanges.  Format:  12 


Record: 

(CHARLGj 

Character i st ic  length  of  dispersive 
exchange  pairings.  Format:  12 

Record : 

( DSP  G  j 

Eddy  dispersion  coefficients  for 
dispersive  exchange  pairings. 

Format  12 

Record : 

iFROCGj 

Organic  carbon  content  of  compartment 
sediments.  Format:  5 

Record : 

iCECG) 

Cation  exchange  capacity  of  sediments 
in  each  compartment.  Format:  5 

Record: 

lAECGj 

Anion  exchange  capacity  of  sediments  in 
each  compartment.  Format:  5 

Record: 

(TCELG) 

Average  temperature  of  ecosystem 
compartments.  Format:  5 

Record : 

( PHG) 

Negative  value  of  log  of  temporal 
average  of  [H+]  concentration  for  each 
compartment.  Format:  5 

Record : 

(POHG) 

Negative  value  of  log  of  temporally 
averaged  lOH-]  concentration  for  each 
compartment.  Format:  5 

Record: 

v0XRAQGj 

Molar  concentration  of  environmental 
oxidants  in  each  ecosystem  compart¬ 
ment  Format:  5 

Record : 

(BIOMSGj 

Total  actively  sorbing  biomass  in  each 
compartment.  Format:  5 

Record : 

(PLRAG) 

Planktonic  fraction  of  total  biomass  in 
each  compartment.  Format:  5 

Record : 

(BIOTMG) 

Biotemperature  in  each  compartment. 
Format:  5 

Record: 

(BACTOGj 

Bacterial  population  density  in  each 
compartment.  Format:  5 

Record: 

(ACBAXG) 

Proportion  of  bacterial  population  that 
actively  degrades  toxicant.  Format:  5 

Record: 

(K02G) 

Rearation  parameter  in  each  compart¬ 
ment.  Format:  5 

Record:  (CMPETGj  Single-valued  zenith  light  extinction 
coefficient  for  water  columns,  dummy 
variable  for  benthic  compartments. 
Format:  5 

Record:  (DQCGj  Dissolved  organic  carbon  concentration 
in  water  column  compartments,  dummy 
variable  in  benthic  compartments. 
Format:  5 

Record:  (CHLG)  Concentration  of  chlorophyll  and 

chlorophyll-like  pigments  in  water 
column  compartments,  dummy  variable 
in  benthic  compartments.  Format:  5 
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